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A  COMPUTER  BASED  APPROACH  TO  ESTIMATING  RUNOFF  CURVE 


NUMBERS  USING  LANDSAT  DATA 
T.  R.  Bondelid,  T.  J.  Jackson  amd  R.  H.  McCuen 

ABSTRACT 

Several  investigations  have  shown  the  potential  of  using  remotely  sensed 
data  in  hydrologic  modeling.  One  of  the  most  promising  is  in  the  estimation 
of  the  land  cover  in  the  computation  of  the  Soil  Conservation  Service  Runoff 
Curve  Number.  Each  of  these  studies  have  shown  these  data  to  be 
cost-effective . 

This  report  describes  a  computer  based  procedure  for  estimating  watershed 
Curve  Numbers  using  remotely  sensed  data.  It  is  a  linkage  of  some  previously 
developed  package  programs  and  hew  procedures  that  have  been  streamlined  for 
this  particular  application, 
types  of  data  could  be  used. 


Landsat  data  are  emphasized,  however,  other 
Examples  are  presented  on  each  aspect. 
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INTRODUCTION 

Recent  research  has  indicated  that  Soil  Conservation  Service  (SCS)  runoff 
Curve  Numbers  (CN)  can  be  accurately  estimated  using  Landsat  imagery 
(Bondelid,  et  al.,  1980;  Ragan  and  Jackson,  1980;  Slack  and  Welch,  1979).  The 
estimation  of  CN's  using  Landsat  are  both  time  and  cost  effective  because  the 
Landsat  data  are  in  machine  readable  form.  Therefore,  digital  computers  can 
be  used  for  many  of  the  required  computations.  An  efficient  CN  estimation 
procedure  has  been  developed  in  the  course  of  a  study  of  three  watersheds  in 
southeastern  Pennsylvania  (Bondelid,  et  al ,  1980).  The  procedure  is  designed 
to  meet  four  principle  criteria:  (1)  the  computations  should  be  compatable  for 
computer  solution;  (2)  the  procedure  should  be  straightforward  and  not  require 
a  high  level  of  expertise  in  either  remote  sensing  or  computer  programming; 

(3)  the  procedure  should  be  versatile  so  that  it  can  be  modified  and  expanded 
for  use  in  a  variety  of  situations;  and,  (4)  the  procedure  should  be 
accessible  to  SCS  hydrologists  and  other  scientists  and  not  require  special 
equipment  that  is  not  normally  available  to  them.  The  objective  of  this 
report  is  to  present  an  efficient  procedure  that  satisfies  these  four  criteria. 

The  procedure  has  been  tested  and  is  fully  operational  on  the  United 
States  Department  of  Agriculture  (USDA)  IBM  370  computer  system  operated  by 
the  Washington  Computer  Center  (WCC).  The  Landsat  analyses  are  performed 
using  the  Office  of  Remote  Sensing  for  Earth  Resources  (ORSER)  of  Pennsylvania 
State  University  computer  software  package  described  by  Borden,  et  al  (1975). 
The  CN  analyses  are  performed  using  a  computer  program,  called  CNPROG,  that  is 
designed  for  performing  efficient  CN  evaluations  using  a  digital  land  cover 
data  base,  such  as  can  be  produced  using  Landsat.  A  full  description  of  the 
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capabilities,  input  requirements,  and  output  products  of  CNPROG  is  included  in 
this  report.  Examples  of  all  of  the  types  of  ORSER  and  CNPROG  runs  are  also 
included.  The  ORSER  users  manual  (Borden,  et  al,  1975)  contains  full 
descriptions  of  these  procedures. 

COMPUTERIZED  CURVE  NUMBER  ESTIMATION 

Runoff  curve  numbers  are  a  function  of  land  cover,  soil  type,  and 
antecedent  soil  moisture  conditions.  Soil  types  are  classified  into  four 
hydrologic  soil  groups  according  to  the  infiltration  capacity  of  the  soil. 

The  four  soil  groups  are  labeled  A,  B,  C,  and  D;  group  A  soils  have  the 
highest  infiltration  rates  while  group  D  soils  have  the  lowest  rates. 

Soils  data  are  generally  spatial  in  nature,  in  that  a  particular  soil  type 
will  occupy  a  particular  region  or  regions  within  the  watershed  under  study. 
Also,  for  hydrologic  modeling  purposes,  the  CN's  must  be  estimated  for  each 
individual  watershed.  The  three  spatial  parameters  of  land  cover,  soil 
groups,  and  subwatersheds  must  be  coordinated  in  producing  the  weighted 
average  CN  for  a  watershed. 

The  basic  function  of  CNPROG  is  to  transform  spatial  soils  and 
subwatersheds  boundaries  to  a  digital  data  base  that  has  a  point  by  point 
correspondence  to  the  digital  land  cover  data  file.  The  program  can  then 
evaluate  runoff  CN's  on  a  point  by  point  basis  and  compute  the  CN's  on  a 
subwatershed  basis.  The  program  also  has  the  capability  of  generating 
subwatershed  and  soils  maps.  These  maps  are  useful  for  checking  the  accuracy 
of  the  input  soil  and  subarea  boundary  point  data. 

Figure  1  presents  a  flow  chart  of  the  computerized  CN  estimation 
procedure.  The  three  basic  steps  in  the  procedure  are:  (1)  development  of 
watershed  data  files;  (2)  development  of  Landsat  land  cover  data  files;  and, 

(3)  integration  of  watershed  and  land  cover  data  to  produce  the  CN  estimates. 
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The  basic  data  requirements  for  performing  the  procedure  are: 

.  Landsat  computer  tapes  and  imagery  of  the  region 

.  United  States  Geological  Survey  (USGS)  maps  of  the  region  at  both 
1:25,000  and  1:1,000,000  scales 

.  A  map  of  the  watershed  and  subwatershed  boundaries  at  a  1:24,000  scale 

.  Either  a  map  of  the  hydrologic  soils  groups  or  a  table  that  gives  the 
percent  of  each  soil  group  for  each  subwatershed 

The  remaining  sections  of  this  report  will  essentially  follow  the  sequence 
shown  in  the  flow  chart  of  Fig.  1.  The  general  procedure  can  be  executed  on 
any  computer  system  that  can  support  a  Landsat  analysis  package  and  CNPROG. 

.  .  .  if'  i 

The  examples  m  this  report  will  use  the  ORSER  package  and  CNPROG  on  the  USDA 
IBM  370. 

The  left  hand  side  of  Fig.  1  shows  the  general  sequence  of  steps  necessary 
for  developing  the  subwatershed  boundary  point  data  and  the  soil  group  data. 
CNPROG  is  used  to  transform  boundary  point  data  to  a  digital  data  base  that  is 
compatible  with  the  Landsat  land  cover  data  base.  Subwatershed  and  soil  group 
maps  can  be  produced  by  CNPROG  for  checking  the  accuracy  of  the  boundary  point 
data.  Finally,  CNPROG  is  used  for  integrating  the  subwatershed,  soil  group, 
and  land  cover  data  bases  co  compute  CN's. 

The  remaining  sections  of  this  report  present  the  basic  principles  and 
procedures  necessary  for  performing  the  sequences  shown  in  Fig.  1  including 
Landsat  analysis  procedure  with  examples  of  each  type  of  computer  run  that 
will  be  needed,  the  steps  necessary  for  generating  the  subwatershed  and  soil 
group  input  data  files,  a  users  manual  for  CNPROG,  how  to  process  the 
subwatershed  and  soil  group  input  data  files  and  integrate  these  data  files 
with  the  land  cover  files  for  computing  CN's,  and,  finally,  a  complete  example 
CN  analysis  using  CNPROG. 
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FIGURE  1.  FLOW  CHART  OF  LANDSAT  CURVE  NUMBER  ANALYSIS  PROCEDURE 
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OVERVIEW  OF  LAND SAT  IN  ANALYSIS 

As  shown  in  Figure  1,  there  are  basically  two  parallel  processes  necessary 
for  using  Landsat  to  estimate  CNs.  The  ORSER  computer  package  is  used  for  the 
Landsat  processing  and  CNPROG  is  used  for  the  watershed  data  processing. 

CNPROG  is  also  used  for  the  integration  of  the  Landsat  and  watershed  data 
files  . 

Computer  Requirements 

Because  the  procedure  is  oriented  primarily  to  the  needs  of  SCS 
hydrologists,  the  example  input  files  include  the  necessary  IBM  Job  Control 
Language  (JCL)  statements  for  use  on  the  USDA  IBM  370  computer  system.  An 
understanding  of  IBM  JCL  is  not  necessary  to  understand  the  Landsat  CN 
analysis  procedure.  A  simple  rule  for  those  users  unfamiliar  with  IBM  JCL  is 
that  all  lines  that  begin  with  a  slash  are  JCL  statements;  all  other  lines  are 
actual  input  data.  A  programming  consultant  may  need  to  be  contacted  if  the 
user  is  unfamiliar  with  IBM  JCL. 

The  Landsat  CN  analysis  procedure  is  designed  to  be  easily  adapted  to 
computer  systems  other  than  IBM.  The  basic  requirements  are  that  a  Landsat 
analysis  package  that  produces  line  printer  maps  and  CNPROG  be  accessible  to 
the  user.  CNPROG  is  written  in  FORTRAN  IV  so  that  it  is  essentially  machine 
independent.  The  Landsat  analysis  package  will  require  minimal  modifications; 
these  modifications  are  explained  in  the  CNPROG  User's  Guide,  which  is 
contained  in  this  manual. 

Using  the  ORSER  Package 

The  ORSER  program  package  is  a  comprehensive  set  of  programs  for  analysis 
of  digital  remote  sensing  data.  Although  it  is  not  necessary,  it  is 
recommended  that  the  user  obtain  access  to  a  copy  of  the  ORSER  User's  Manual 
for  the  Landsat  analysis. 
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All  of  the  ORSER  computer  runs  used  in  the  CN  analysis  are  straightforward 
and  require  relatively  few  input  cards  to  perform  the  Landsat  phase  of  the 
procedure.  The  example  input  files  contained  in  this  manual  can  be  used  with 
only  the  obvious  changes  in  various  numerical  values  and  data  file  names. 

The  ORSER  input  cards  consist  of  English  language  commands  and,  if 
necessary,  input  data  values.  For  instance,  the  command  CORRECT  specifies 
that  the  Landsat  data  will  be  geometrically  corrected  during  program 
execution.  Generally,  only  the  commands  contained  in  the  example  input  files 
will  be  necessary  for  the  Landsat  analysis 
Using  CNPROG 

CNPROG  is  a  user  oriented  package  that  can  be  used  in  either  an 
interactive  or  a  batch  mode.  CNPROG  allows  the  user  to  specify  the  logical 
unit  numbers  and  formats  of  the  required  data  files  during  execution.  The 
data  files  that  are  required  for  a  particular  execution  depend  upon  the  CNPROG 
options  used.  The  data  file  requirements  will  be  discussed  in  subsequent 
sections . 

The  basic  input  sequence  consists  of  command  words  that  are  followed  by 
input  data.  For  instance,  the  command  INFILE  is  used  to  specify  an  input  data 
file  unit.  After  the  command  INFILE  is  typed,  the  program  prompts  the  user 
with  a  message  stating  that  the  input  file  unit  should  be  entered  in  columns 
1-2.  The  value  entered  is  then  read  and  checked  for  obvious  input  errors.  If 
the  input  value  is  either  nonnumer ical  or  less  than  or  equal  to  zero,  a 
message  is  printed  that  indicates  an  error.  Then  the  program  again  prompts 
the  user  for  the  input  file  unit  number.  If  the  input  value  is  accepted,  then 
the  program  will  prompt  the  user  for  a  new  command.  This  process  is  continued 
until  all  desired  processing  is  completed. 
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In  an  interactive  mode,  the  program  prompts  and  error  messages  guide  the 
user  through  the  program  execution.  Also,  an  incorrect  data  entry  will 
generally  not  cause  an  abortion  of  the  entire  program  execution.  In  a  batch 
mode,  the  prompts  and  error  messages  are  suppressed.  If  an  input  error  is 
detected,  then  program  execution  continues  until  actual  data  file  processing 
starts.  At  that  point,  if  errors  have  been  detected,  error  messages  are 
printed  and  data  file  processing  is  not  performed. 

CNPROG  Processors 

CNPROG  consists  of  basically  three  processor  options.  The  first  processor 
is  the  file  generator  and  is  invoked  by  the  command  FILEGEN.  The  file 
generator  option  is  used  for  producing  soils  and  subarea  files  that  have  a 
point  by  point  correspondence  to  the  land  cover  data  files.  These  files  can 
be  either  saved  for  use  in  subsequent  executions  or  used  in  the  same 
execution.  The  second  processor  is  the  map  generator  that  is  invoked  by  the 
command  MAP GEN.  This  processor  is  used  for  producing  subarea  and  soils  maps 
from  files  produced  by  FILEGEN.  The  third  processor  is  the  CN  calculator, 
which  is  invoked  by  the  command  CNCALC.  This  processor  utilizes  the  subarea 
and  soils  files  generated  by  FILEGEN  and  the  land  cover  data  files  generated 
by  the  land  cover  classifier.  The  CNCALC  processor  computes  areas,  CN's,  land 
cover  percentages,  and  soil  type  percentages  for  each  subarea  and  for  the 
watershed  as  a  whole. 

All  of  the  processors  can  be  used  in  one  execution  or,  for  instance,  only 
the  FILEGEN  processor  can  be  used  in  one  program  execution.  The  generated 
files  in  the  first  execution  can  be  used  in  subsequent  executions  in  which  the 
MAPGEN  and/or  CNCALC  processors  are  used. 
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FILEGEN  Processor 

The  file  generation  processor  (FILEGEN)  converts  the  subarea  and  soils 
boundary  point  data  to  files  that  have  a  point  by  point  correspondence  to  the 
land  cover  data  files  produced  by  the  land  cover  classifier.  Sometimes  the 
map  image  for  an  area  being  studied  is  wider  than  the  line  printer  paper.  In 
these  cases,  the  map  Image  consists  of  two  or  more  map  strips,  which  can  be 
placed  side  by  side  to  produce  the  entire  map.  The  land  cover  map  image  file 
thus  is  either  sequenced  by  these  map  strips  or  there  can  be  a  separate  file 
for  each  map  strip.  The  files  generated  by  FILEGEN  must  also  be  sequenced  by 
map  strips  corresponding  to  the  land  cover  map  strips. 

The  map  images  produced  by  land  classifiers  can  be  thought  of  as  a 
coordinate  system  in  which  each  print  line,  referred  to  as  a  scan  line,  is 
numbered  sequentially  from  top  to  bottom.  Each  column,  referred  to  as  an 
element,  is  numbered  sequentially  from  left  to  right.  Each  pixel  on  the  map 
is  thus  defined  by  scan  line  and  element  values,  which  are  produced  by  the 
land  cover  classifier.  Each  land  cover  type  is  assigned  an  alphameric 
character  so  each  scan  line  consists  of  a  string  of  characters  in  which  each 
character  specifies  the  land  cover  for  a  pixel. 

The  land  cover  data  can  be  in  a  sequential  file  that  is  (n  x  m)  records 
long,  in  which  n  is  the  number  of  scan  lines  and  m  is  the  number  of  map 
strips.  The  first  n  records  contain  the  land  cover  data  for  the  first  map 
strip,  the  next  n  records  are  for  the  second  map  strip,  etc.  If  each  map 
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strip  is  on  a  separate  file,  then  there  are  m  files  with  each  file  being  n 
records  long.  The  number  of  elements  in  each  record  is  equal  to  the  width  of 
the  particular  map  strip. 

The  files  produced  by  FILEGEN  are  analogous  to  the  land  cover  map  files; 
they  are  organized  sequentially  by  scan  line  and  map  strips.  The  user 
generally  does  not  need  to  understand  this  file  structure  because  CNPROG 
internally  handles  the  file  reading  and  writing. 

MAPGEN  Processor 

The  MAPGEN  processor  generates  line  printer  map  images  from  the  FILEGEN 
scan  line  files.  These  maps  use  the  same  line  and  element  coordinate  system 
that  is  used  for  the  land  cover  file.  The  maps  are  useful  for  checking  the 
accuracy  of  the  subarea  and  soils  boundary  data.  CNPROG  automatically  assigns 
a  unique  map  symbol  to  each  subarea. 

CNCALC  Processor 

The  CN  calculations  are  performed  by  the  CNCALC  processor.  CNCALC  uses 
the  land  cover  file  and  FILEGEN  scan  line  files  for  computing  the  CN's.  The 
CN's  are  calculated  for  each  point  for  each  subarea  and  takes  a  weighted 
average  for  computing  the  overall  CN  for  each  subarea. 

The  CN  for  a  particular  pixel  is  based  on  an  input  CN  table.  The  table 
contains  the  CN  values  for  the  four  hydrologic  soil  groups  for  each  alphameric 
land  use  symbol.  If  the  land  cover  for  a  point  is  unclassified  (i.e.,  the 
symbol  for  that  pixel  is  a  blank  space),  then  the  pixel  is  not  included  in  the 
CN  computation.  The  pixel  is  included  in  the  determination  of  total  area. 


11 


LANDSAT  LAND  COVER  CLASSIFICATION 

The  basic  objective  of  the  Landsat  analysis  is  to  produce  a  land  cover  map 
that  can  be  used  for  estimating  CN's.  The  land  cover  classifications  must  be 
sufficiently  general  to  be  identifiable  with  Landsat  while  also  containing 
enough  detail  to  adequately  define  the  CN  values-  Table  1  contains  an  example 
of  a  set  of  Landsat  land  cover  categories  and  the  corresponding  CN  values  for 
the  four  hydrologic  soil  groups.  This  table  was  used  successfully  in  a  study 
in  southeastern  Pennsylvania  (Bondelid,  et  al . ,  1980).  Table  1  is 
presented  for  illustrative  purposes  only;  the  land  cover  classifications 

TABLE  1.  Example  Landsat  Curve  Number  Table 


Land  Cover 

A 

Soil 

B 

Group 

C 

D 

Woods 

25 

55 

70 

77 

Agriculture 

64 

75 

83 

87 

Residential 

60 

74 

83 

87 

Highly  Impervious 

SO 

93 

94 

95 

Water 

98 

98 

98 

98 

used  in  a  particular  region  may  vary  depending  on  the  dominate  land  cover 
types  in  the  region. 

An  understanding  of  the  basic  properties  of  Landsat  data  is  essential  for 
understanding  the  Landsat  land  cover  classification  procedure.  The  Landsat 
data  is  available  in  two  forms?  (1)  Landsat  imagery,  in  which  the  data  are 
presented  as  color  or  black  and  white  pictures  at  a  scale  of  1:1,000,000}  or, 
(2)  the  Landsat  digital  data,  which  are  contained  on  computer  compatible 
magnetic  tapes. 
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Landsat  Data 

The  basic  Landsat  classification  procedure  is  a  four  step  process: 
procuring  the  Landsat  data,  preliminary  data  processing,  determining  the  land 
cover  classifications,  and  producing  a  land  cover  map.  The  final  land  cover 
map  is  produced  as  a  line  printer  map  and  also  in  a  machine-readable  form  for 

computing  CN's  using  CNPROG. 

The  Landsat  satellites  circle  the  globe  in  a  circular,  near-polar, 
sun-synchronous  orbit  so  that  the  same  point  on  the  Earth's  surface  is  viewed 
every  18  days  at  the  same  time  of  day.  The  sensing  system  that  is  used  for 
the  CN  analyses  is  the  multispectral  scanner  (MSS),  which  acquires  images  in 
four  spectral  bands  of  the  electromagnetic  spectrum.  These  four  bands  ar?: 

1.,  Band  4,  the  green  band,  between  0.5  and  0.6  micrometers  (m). 

This  band  emphasizes  sediment-laden  and  shallow  waters. 

2.  Band  5,  the  red  band,  between  0.6  and  0.7  m.  This  band 
emphasizes  man-made  features  such  as  ufban  and  rural  settlements. 

3.  Band  b,  a  near  infrared  ( IR)  band,  between  0.7  and  0.8  m. 

This  band  emphasizes  vegetation  and  landforms. 

4.  Band  7,  another  near  IR  band,  between  0.8  and  1.1  m.  This 
band  also  emphasizes  vegetation  and  landforms,  and  provides  the 
best  haze  penetration. 

The  Landsat  data  are  separated  into  sections  called  "scenes".  Each  scene 
is  approximately  185  km.  wide  by  185  km.  long.  The  Landsat  computer 
compatible  tapes  (CCT)  are  stripped  into  four  separate  files,  with  each  strip 
being  46.25  km  wide  by  185  km.  long. 

Each  digital  picture  element  (pixel)  is  57  m  wide  by  79  m  long.  The 
Landsat  CCT's  contain  the  spectral  response  values  in  the  four  bands  for  each 
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4  & 

pixel  in  the  scene.  The  components  of  Landsat  scenes  are  illustrated  in 
Figure  2.  The  rows  of  pixels  are  referred  to  as  "scan  lines"  and  the  columns 
are  referred  to  as  elements.  The  lines  are  numbered  consecutively  from  north 
to  south  and  the  elements  are  numbered  from  west  to  east. 

Landsat  classification  of  land  cover  is  based  on  the  premise  that  each 
land  cover  type  ha3  a  specific  range  of  spectral  responses  in  one  or  more 
channels.  For  instance,  because  water  has  very  low  IR  reflectance,  all  pixels 
with  low  reflectance  values  in  bands  6  or  7  can  be  classified  as  water.  The 
appropriate  spectral  response  ranges  must  be  determined  for  each  scene  being 
analyzed  because  of  variables  such  as  sun  angle  and  atmospheric  haze. 

Information  on  ordering  Landsat  data  can  be  obtained  from  the  Earth 
Resources  Observation  Systems  (EROS)  Program  center  administered  by  the  USGS. 

A  packet  of  ordering  information  can  be  obtained  from: 

Customer  Relations 
EROS  Data  Center 
Sioux  Falls,  SD  57198 
(605)  594-6511  ext.  151 

'  •% 

Preliminary  Processing 

After  the  Landsat  data  have  been  obtained,  some  preliminary  processing  is 
necessary.  The  preliminary  processing  includes:  (1)  subsetting  the  study  area 
from  the  full  Landsat  scene;  (2)  geometrically  correcting,  rotating,  and 
rescaling  the  data;  and,  (3)  producing  a  "brightness"  map  of  the  study  area. 
The  ORSER  programs  that  are  used  for  the  preliminary  processing  are  SPAN, 
SUBGM,  and  NMAP. 
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FIGURE  2.  LANDS AT  SCENE  COMPONENTS 
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Subsetting  the  Study  Area 

The  area  under  study  for  CN  analysis  will  generally  be  much  smaller  than 
the  full  Landsat  scene.  For  efficient  processing,  a  subset  tape  that  contains 
data  for  only  the  general  study  area  should  be  created.  The  Landsat  data  can 
also  be  geometrically  corrected,  rescaled  to  1:24,000,  and  rotated  to  true 
north  during  the  subsetting  procedure. 

The  area  to  subset  is  defined  by  a  rectangular  block  in  which  the  line  and 
element  coordinates  of  the  block  are  input  to  the  appropriate  ORSER  program. 
The  coordinates  of  the  block  can  be  determined  by  overlaying  a  transparent 
1:1,000,000  scale  map  over  the  Landsat  imagery.  The  rectangular  block  for  the 
study  area  can  then  be  outlined  on  the  Landsat  image.  The  coordinates  (line 
and  element  values)  of  the  block  can  be  easily  determined  because  the  number 
of  lines  and  elements  in  the  entire  scene  is  known.  For  instance,  if  the  left 
side  of  the  block  is  located  35%  of  the  total  width  away  from  the  left  edge  of 
the  scene,  then  the  approximate  element  value  of  the  left  side  of  the  block  is 
0.35  times  the  total  width  in  pixels  (approximately  3,230),  which  is  equal  to 
1130. 

If  the  block  to  be  subset  overlaps  the  border  between  two  adjacent  strips, 
then  a  temporary  subset  tape  must  be  created  using  the  ORSER  program  SPAN. 

SPAN  will  produce  a  single  ORSER  subset  tape  using  the  two  Landsat  files  for 
the  adjacent  strip.  This  temporary  tape  is  then  used  in  generating  the  final 
subset  tape.  If  the  block  is  entirely  on  one  strip,  then  SPAN  does  not  need 
to  be  run. 

Correcting  and  Rescaling 

The  final  subset  tape  can  be  produced  using  the  ORSER  program  SUBGM.  This 
program  can  subset,  geometrically  correct,  rotate,  and  rescale  the  data.  The 
data  must  be  geometrically  corrected  because  the  Landsat  image  is  skewed  due 
to  the  rotation  of  the  Earth.  The  data  must  also  be  rotated  to  the  true  north 
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because  the  Landsat  path  is  not  exactly  north-south.  The  data  must  also  be 
rescaled  so  that  the  line  printer  maps  will  be  at  a  scale  of  1:24,000. 
Producing  a  "Brightness"  Map 

After  the  subset  tape  has  been  generated,  a  map  that  shows  the  overall 
patterns  of  "light"  and  "dark"  reflecting  areas  should  be  produced.  Major 
landforms  such  as  mountain  ranges  and  water  bodies  will  appear  as  dark  areas 
on  the  map.  These  landforms  can  then  be  matched  up  with  the  1:24,000  USGS 
topo  maps  to  get  accurate  alignment  betwen  the  ORSER  map  and  the  topo  maps. 
Example  Input  Files 

Figure  3  shows  an  example  input  file  for  running  the  SPAN  program. 

Fortran  units  8  and  9  contain  the  two  strips  to  be  spanned  and  unit  10  is  for 
the  new  tape  which  contains  the  spanned  block.  The  "SPAN"  input  card 
specifies  that  the  block  written  to  unit  10  will  contain  lines  1030  to  1430 
and  elements  1580  to  2025.  The  "CHANNELS"  card  specifies  that  all  four 
channels,  that  is,  all  four  Landsat  bands,  will  be  included  on  the  new  tape. 
The  "END"  card  specifies  the  end  of  input  data. 

Figure  4  shows  an  example  SUBGM  input  file.  The  new  tape  produced  by  the 
example  run  of  Figure  3  is  used  as  the  input  tape.  The  "BLOCK"  card  specifies 
that  a  block  containing  lines  1030  to  1430  and  elements  1580  to  2025  will  be 
processed.  The  "CORRECT",  "ROTATE",  "CUBIC",  and  "ALLOW"  cards  specify  that 
the  new  tape  will  contain  geometrically  corrected,  rotated,  and  rescaled 
Landsat  data. 

Figure  5  shows  an  example  NMAP  run  using  the  tape  created  by  running 
Figure  4.  Note  that  units  6  and  9  are  used  for  line  printer  output.  Each 
unit  prints  one  map  strip;  no  more  than  two  map  strips  can  be  printed  in  a 
given  run  of  NMAP. 


//SEANHTB1  JOB  <♦♦♦♦♦♦♦♦♦♦, RJ0£9>  »  ' BDNBEL IB-'* MSSLEV'EL=  <1» 1>  * 

//  CL ASS=B  >  T I ME= < *  2  0>  > PRTY=£ 

/♦ROUTE  PRINT  RMT£9 

//JOEL  I B  HD  DSN=SEANH  02 .  DRSER£ » B I SP-SHR 
//STEP  1  EXEC  PGM=SPAN  ?  REG  I  DN=£  0  OK 
//FT06F001  BB  SYSOUT=fl 

//FT0SF001  BB  UNIT=TRPE9?  V0L=SER=A1£345?  BSN=SERNHO£.  TAPE!  .  BATA.* 

//  DISP=  <OLB j  KEEP) j  LABEL=  <£*  SL? >  HO 

//FT09F001  BB  UNIT=TAPE9> V'DL=SER=B1£345.* BSN=SEANHO£. TAPER. BRTA> 

//  LRBEL= <1 , SL>  > 1N>  ? BISP= <DLI»?  KEEP> 

//FT1 OF 001  BB  UNXT=TAPE9? BISP= <NEW*  KEEP* BELETE> ?  BSN=SERNHO£. TAPES. BATf 
//  BCB= <RECFM=VBS?  LRECL=3696>  BLKSI2E=3700> .*  LABEL=  <1 , SL> 

//FT05F001  BB  ♦ 

SPAN  1030  1430  1580  £0£5  1  1 

CHANNELS  1  £  3,4. 

END 

/♦ 


.  FIGURE  3.  EXAMPLE  SPAN  INPUT  FILE 


✓✓SEflNHTBi  JDB  <♦♦♦♦♦♦♦♦♦♦,  R.J02SO ? 'BDNBELIB' ? MS6LEVEL= <1 * 1> > 

✓✓.  CLflSS=M T I  ME=  1  0  >  PRT Y=  1 
✓♦ROUTE  PRINT  RMT29 

✓✓JOEL I B  BB  BSN=SE ANH  02 . 0RSER2 .< B I SP=SHR 
✓/STEP 1  EXEC  PGM=SUBGM .« REG I DN=2 0 OK 
✓✓FT06F001  BB  SYSQUT=fi 

✓✓FT  OSF  0  0 1  BB  UN  I  T=TAPE9  >  VOL=SER=C  1  £345 BSN=SERNH  02 .  TRPE3 .  BflTfi  ? 

✓✓  BISP= <DLB>  KEEP> f LREEL=  Cl » SL» > IN> 

/✓FT09F001  BB  UNIT=TRPE9*  BISP=  <NEW»  KEEP?  BELETE)  >  BSN=SERNH 02 .  TRPE4.  DATA 
✓✓  BCE= <RECFM=VBS * LRECL=3696 > ELKS I ZE=37  0 0>  ?  LREEL= < 1 » SL> 

✓✓FT05F001  BB  ♦ 

BLOCK  •  1030  1430  1580  2025  1  1 

CHANNELS  1234 
CORRECT 
ROTATE 
CUBIC 
ALLOW 
ENB 
✓♦ 


FIGURE  4 


EXAMPLE  SUBGM  INPUT  FILE 


//SEfiNHTBl  JOB  <♦♦*♦♦♦*♦♦♦, RJ0£9> » 'EDNDELID'* MSGLEVEL= <1 , 1) , 

//  CLASS=B ,  1 1  ME=  < £  0>  >  PRTY=2  '  •  . 

/♦ROUTE  PRINT  RMT£9 

//JOEL  I E  HD  DSN=SERNH  OS .  .ORSERS  ?  D I  SP=SHR 
//STEP!  EXEC  F‘GM=NMRP .<  REG  I  ON=£ 0  OK 
//FT06F001  DD  SYSDUT=R 

//FT08F001  DD  UNIT=TRPE9> VOL=SER=D 12345,  DSN=SEANH02. TRFE4. DATA, 

//  DISF  -  COLD.-  KEEP)  i-  LREEL=  <1 >  SL, ?  I  TO 

//FT09F001  DD  SYSDUT=fl» DCB= <RECFM=Ffl, ELKSIZE=133> 

//FT05F001  DD  ♦ 

BLOCK  1030  1430  1580  1780  1  1 

CHANNELS  1  £  3  4 

END 

/♦ 


FIGURE  5 


EXAMPLE  NMAP  INPUT  FILE 
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Land  Cover  Classification 

The  spectral  responses  for  the  desired  land  cover  classifications  must  be 
determined  for  each  Landsat  scene  being  analyzed.  The  first  step  is  to  select 
"training"  sites  in  which  a  particular  land  cover  type  is  dominate.  An  ideal 
training  site  contains  only  one  particular  land  cover  type.  The  spectral 
responses  of  the  training  sites  can  be  used  to  generate  the  ranges  of  spectral 
responses  applicable  to  each  land  cover  type.  These  ranges  of  spectral 
responses  can  then  be  used  to  define  the  land  cover  classifications  for 
producing  a  land  cover  map. 

The  statistical  analysis  of  the  training  sites  can  be  performed  using  the 
modified  ORSER  program  STATS2.  This  program  is  a  modified  version  of  the 
STATS  program  in  which  the  capability  of  producing  two  dimensional  histograms 
of  the  spectral  responses  of  two  channels  has  been  added. 

The  two  dimensional  histogram  option  in  STATS2  is  invoked  by  using  the 
2DHIST  control  card.  2DHIST  is  specific  to  the  STATS2  program  on  the  USDA 
computer  and  does  not  exist  in  the  standard  ORSER  package.  The  format  of  the 
2DHIST  control  card  is  presented  in  Table  2.  An  example  STATS2  input  file  is 
shown  in  Figure  6.  The  "SIGN"  card  specifies  that  a  summary  of  the  spectral 
signatures  for  each  training  area  will  be  output  to  unit  11.  The  "HIST"  card 
specifies  that  one -dimensional  histograms  for  channels  2  and  4  will  be 
produced.  Each  set  of  "CATEGORY",  "AREA",  and  "ENDAREA"  cards  specifies  a 
training  site.  The  "CATEGORY"  card  specifies  the  training  site  index  number 
and  name.  Each  "AREA"  card  specifies  the  line  and  element  values  of  a  point 
on  the  training  site.  The  "AREA"  cards  must  be  input  sequentially  clockwise 
or  counter-clockwise  around  the  training  site  polygon.  The  "ENDAREA"  card 
specifies  the  end  of  input  polygon  points  for  the  training  site. 
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TABLE  2.  2DHIST  Control  Land  Format 

Cols . 

Description  of  Contents 

1-6 

2DHIST 

10 

Flag  for  plot  size: 

1  =  60  characters  wide  by 

40  characters  high 

2  =  100  characters  wide  by 

50  characters  high 

11-26 

8  two-column  integer  fields  specifying  the  channels  to 

be  used  for  the  histograms.  Each  histogram  requires 

two  fields;  the  first  field  specifies  the  channel  to 

be  used  for  the  abscissa  and  the  second  field 

specifies  the  channel  to  be  used  for  the  ordinate.  Up 

to  four  different  histograms  can  be  specified. 

31-50 

4  five-column  real  fields  which  specify  the 

approximate  ranges  of  values  to  be  displayed.  The 

values  are  input  in  the  order:  minimum  x  value, 

maximum  x,  minimum  y,  maximum  y.  If  these  values  are 

left  blank,  then  the  program  will  determine 

appropriate  values  for  each  histogram.  Specifying  the 

x  and  y  ranges  is  useful  when  comparing  histograms  for 

different  training  sites  because  the  scales  will  be  the 

same  for  all  two  dimensional  histograms. 

V/SERNHTB6  JOB  <4071 090299* RJ029>  *  'BOMBELID' *MSSLEVEL“<1 *  i>  * 

-V/  CLASSES *TIME=<*20> 

/*ROUTE  PRINT  RMT29 

//JOBLIB  BB  DSN—SEANH02  .  DRSER2 *  B 1 SP=SBR 
//STEP!  EXEC  P6M=STRTS2 * REGI DN=3 0 ©K 
//FT06F00!  BB  SYSDUT=A 

//FT09P00S  BD  UNIT-TAPE9*VDL==SER==028758*BSN=ARS38.PAOS1R^ 

✓✓  LflBEL= <1*SL*?IN>*BI SP=  <DLB *  K£EP> 

//FT11F001  DB  UNIT«SYSDfiyDSN=SEfiNH02.TRB. STOUT 1 -BATB>*BISP-<NEW*CBTLG> * 
//  BCB-  <RECFM”FB  *  LRECL=8  0  *  BLKS I  ZE=3 1 20>  ?  SPACE^  <3 1  £  0  *  <25  »  5>  *  RLSE> 

//FT05F001  BB 
SIGN  *  11 

CHANNELS  .'2  4-  :  ■  .  . 

HIST  2  4  !  " 

SCHIST  1  2400000©  0.0  59.  O.O  73. 

CATEGORY  1  AGRICULTURE  - 

AREA  906  1798  - 

AREA  906  1808 

AREA  .  911  1808  * 

AREA  911  1798  -  ‘ ^ 

ENDAREA  •  ’  ‘ 

CATEGORY  2  FOREST  1 

AREA  976  1640  .  .  '  ' 

AREA  976  1660  -  ' 

AREA  •  983  1660  '»  * 

AREA  ■  983  1640  '  . .  • 

ENBAREA  ✓  '  .  - 

CATEGORY  3  FOREST  2  _  •  ' 

AREA  942  1860 

BREA  942  1878 

RREA  •  /  952  1878  '  •  • 

AREA  952  I860 

ENBAREA  • 

CATEGORY  4  URBAN  -  - 

AREA  872  1731 

AREA  872  1755 

BREA  880  1755 

AREA  880  1731 

ENBAREA  *' 

CATEGORY  5  RESIB.  MOSTLY 
AREA  890  1730 

AREA  890  1738 

AREA  895  1738 

BREA  895  1730 

ENBAREA 


CATEGORY  6 

AREA 

AREA 

AREA 

AREA 

ENBAREA 


MINE  TAILINGS  POND 
862  1714 
862  1716 
865  1716 
865  1714 


END 


FIGURE  6 


EXAMPLE  STATS2  INPUT  FILE 
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Parallelapiped  classification  method, 

•  &" 

# 

All  methods  of  classifying  Landsat  data  are  ba^o  on  the  premise  that  a 
particular  land  cover  type  has  a  distinct  range  of  spectral  responses.  The 
parallelapiped  classification  (PPD)  method  utilizes  rectangular  ’'boxes''  to 
specify  the  ranges  of  spectral  responses  for  each  land  cover  type.  The  PPD 
method  is  recommended  for  Landsat  CN  analyses  for  four  reasons:  (1)  the  method 
has  been  used  successfully  in  previous  CN  studies;  (2)  the  method  is  easy  to 
understand  and  use;  (3)  the  number  of  unclassified  pixels  can  be  kept  small; 
and,  (4)  the  method  is  relatively  inexpensive  in  terms  of  computer  costs. 

Frequently,  only  two  bands  are  necessary  to  adequately  classify  Landsat 
data;  one  visible  band  and  one  IR  band.  This  is  because  the  spectral 
responses  for  a  particular  land  cover  type  are  very  similar  when  comparing 
either  the  two  visible  or  the  two  IR  bands  to  each  other.  For  instance,  the 
study  in  southeastern  Pennsylvania  used  only  the  data  for  bands  5  and  7  in 
performing  the  PPD  classifications.  The  classification  boundaries  used  in 
that  study  are  shown  in  Figure  7.  Note  that  more  than  one  "box"  was  used  for 
defining  each  land  cover  type. 


LANDSAT  BAND  7  RESPONSE 
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The  ORSER  program  PPD  produces  land  cover  maps  using  the  parallelapiped 
method.  Each  PPD  run  can  produce  a  line  printer  map  that  is  up  to  two  sheets 
wide.  If  the  study  area  is  wider  than  two  sheets,  then  more  than  one  run  of 
PPD  is  required. 

The  two  dimensional  histograms  produced  by  STATS2  are  quite  useful  in 
doing  a  PPD  classification.  The  following  example,  using  data  from  the  study 
in  southeastern  Pennsylvania,  is  presented  to  illustrate  the  PPD 
classification  procedure. 

Figures  8  through  13  show  histograms  for  selected  training  sites.  The 
classification  boundaries  from  Figure  2-5  are  superimposed  on  the  histograms 
to  illustrate  how  the  PPD  limits  can  be  determined  using  the  histograms. 

Figure  6  shows  the  STATS2  input  file  used  for  producing  Figures  7  to  13. 

The  histograms  show  the  frequency  distributions  of  the  spectral  response 
values  in  selected  training  areas.  The  training  areas  used  in  the  examples 
are  agriculture,  forest,  urban,  residential,  and  a  mine  tailings  pond.  Each 
histogram  represents  the  spectral  distribution  for  a  particular  land  cover 
type.  The  PPD  bounds  are  delineated  by  drawing  rectangular  boxes  around  the 
clusters  of  points  on  the  histograms.  Note  that  several  of  the  boxes  drawn  on 
the  example  histograms  are  larger  than  the  clusters  shown.  A  particular  PPD 
"box"  can  frequently  be  enlarged  so  long  as  the  enlargement  does  not  interfere 
with  the  bounds  for  other  land  cover  types. 

The  classification  process  is  often  iterative,  as  shown  in  Figure  1.  The 
land  cover  map  determined  from  the  first  STATS2  analysis  may  not  be  adequate; 
that  is,  there  may  be  regions  of  unclassified  or  mis-classified  land  covers. 

In  such  cases,  additional  training  sites  should  be  analyzed  and  the  PPD  bounds 
adjusted  until  the  resulting  land  cover  map  is  judged  adequate. 
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FIGURE  8 .  EXAMPLE  2-D  HISTOGRAM  FOR  AGRICULTURE 
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FIGURE  10 


EXAMPLE  2“D  HISTOGRAM  FOR  FOREST  2 
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EXAMPLE  2-D  HISTOGRAM  FOR  URBAN'  * 
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EXAMPLE  2-D  HISTOGRAM  FOR  RESIDENTIAL 
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FIGURE  13 


EXAMPLE  2-D  HISTOGRAM  FOR  MINE  TAILINGS  POND 
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Land  Cover  Data  Files 

The  final  land  cover  map  must  be  in  machine-readable  form  for  computing 
the  CN's  with  CNPROG.  The  PPD  program  has  been  modified  to  produce  land  cover 
data  files  that  can  be  used  in  CNPROG,  The  modified  PPD  program  is  called 
PPDSAV.  PPDSAV  writes  the  line  printer  map  images  to  logical  units  20  and  21; 
unit  20  contains  the  second  map  strip.  If  only  one  map  strip  is  produced, 
then  only  unit  20  is  used. 

There  should  be  a  separate  file  for  each  map  strip  in  the  study  area.  The 
control  cards  for  PPDSAV  are  identical  to  the  PPD  control  cards.  An  example 
input  file  for  PPDSAV  is  shown  in  Figure  14.  The  "CATEGORY"  card  specifies 
the  category  (i.e.,  land  cover  type)  number,  map  symbol  for  use  on  the  line 
printer  map,  and  the  category  name.  The  "LOWER"  and  "UPPER"  cards  specify  the 
lower  and  upper  spectral  bounds,  respectively,  for  a  category.  Because  the 
"CHANNELS"  card  specifies  use  of  only  channels  2  and  4,  the  spectral  bounds  on 
each  "LOWER"  and  "UPPER"  card  specifies  the  bounds  for  channels  2  and  4. 

PPDSAV  also  produces  the  line  printer  maps.  The  line  and  element  limits 
of  each  map  strip  should  be  noted  because  this  information  is  required  in 
running  CNPROG.  CNPROG  also  requires  as  input  the  format  specifications  of 
the  land  cover  files.  The  land  cover  file  format  that  should  be  input  to 
CNPROG  when  using  PPDSAV  is  (IX,  15,  2X,  122A1). 

The  generation  of  the  land  cover  data  files  is  the  last  step  of  the 
Landsat  land  cover  analysis.  The  land  cover  data  files  must  be  coordinated 
with  the  watershed  data  files.  The  procedure  for  generating  the  watershed 
data  files  is  explained  later. 
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CNPROG  USER'S  GUIDE 

CNPROG  can  be  used  for  processing  the  watershed  data,  generating  line 
printer  maps,  and  computing  CNs.  The  subwatershed  and  soils  data  file 
requirements  will  be  discussed  in  this  section.  A  full  description  of  the 
remaining  CNPROG  input  requirements,  program  limitations,  and  output  products 
is  also  presented  in  this  section. 

CNPROG  Limitations 

CNPROG  has  upper  limits  for  several  parameters.  These  upper  limits  are 
generally  based  on  limiting  the  computer  storage  requirements,  while  still 
being  able  to  meet  most  of  the  anticipated  requirements.  The  total  map  width 
is  limited  by  the  maximum  number  of  elements  per  strip  and  the  maximum  number 
of  strips.  The  maximum  number  of  elements  allowed  per  map  strip  is  120;  this 
limitation  is  due  to  the  physical  size  of  line  printer  output.  The  maximum 
number  of  map  strips  allowed  is  seven,  so  the  maximum  number  of  elements  (the 
map  width)  is  840.  There  are  no  limits  on  the  number  of  lines  in  a  map. 

There  are  no  limits  on  the  size  of  any  subarea  or  soil  area,  although  the 
maximum  number  of  total  boundary  points  for  all  subareas,  including 
interpolated  points,  is  currently  set  at  6000. 

The  number  of  soil  types  is  limited  to  the  four  hydrologic  soil  types  A, 

B,  C,  and  D.  The  maximum  number  of  subareas  allowed  is  36;  this  limitation  is 
set  by  the  number  of  available  discrete  maps  symbols.  There  are  no  limits  to 
the  number  of  closed  polygons  allowed  per  subarea. 

CNCALC  can  consider  a  maximum  of  ten  land  cover  types.  If  more  than  ten 
land  cover  types  are  needed,  then  the  arrays  in  CNCALC  can  be  easily 
redimensioned  to  whatever  value  is  needed. 
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In  any  program  execution  there  are  no  limits  to  the  number  of  calls  or  the 
sequence  of  calls,  to  any  of  the  three  processors.  Thus,  any  number  of 
analyses,  or  combinations  of  analyses,  can  be  performed  in  one  program 
execution. 

CNPROG  Input 

The  user  controls  the  execution  and  input  sequences  by  entering  simple 
English  language  commands.  CNPROG  is  versatile  in  several  respects.  The 
sequence  of  commands,  and  hence  the  input  sequence,  is,  in  general,  up  to  the 
user.  The  formats  and  logical  unit  numbers  of  input  data  files  are  user 
defined.  Also,  the  program  has  numerous  checks  on  the  validity  and 
consistency  of  the  input  data  so  tl'iat  incorrect  data  entry  will  usually  not 
result  in  a  termination  of  the  entire  run . 

The  CNPROG  commands  and  basic  input  requirements  will  be  detailed  in  this 
section.  The  required  data  files  depend  upon  the  CNPROG  options  selected. 

The  data  files  will  be  discussed  in  the  next  section. 

There  are  two  types  of  commands;  main  commands  and  processor  commands. 

The  main  commands  are  used  to  invoke  the  desired  processor  or  to  stop  program 
execution.  The  processor  commands  are  used  within  a  processor;  each  processor 
has  a  unique  set  of  commands. 

The  basic  input  sequence  is  as  follows:  first,  a  main  command  is  entered 
which  invokes  the  desired  processor.  Next,  a  series  of  processor  commands  and 
data  is  entered  until  the  processing  is  complete.  Control  is  then  returned  to 
the  main  program,  and  either  another  processor  is  invoked,  or  execution  is 
terminated  by  the  command  STOP.  Any  number  of  calls  to  the  processors  can  be 
performed  in  a  single  program  execution. 
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Command  Notation 

Each  command  is  an  alphabetic  string  of  characters  that  starts  in  column 
one.  Several  commands  specify  that  particular  input  data  is  to  be  entered. 

All  integer  input  must  be  right-justif ied ;  the  last  digit  of  the  number  should 
be  in  the  last  column  of  the  data  field  specified.  All  real  valued  input 
should  include  the  decimal  point. 

All  of  the  CNPROG  commands  are  described  below.  The  descriptions  include 
the  purpose  of  the  command,  the  program  response,  and  the  required  user 
response . 

MAIN  Command : 

1.  Command :  FILEGEN 

Purpose :  for  generating  point  record  files  for  subareas  or  soils. 

Program  Response:  calls  the  sequence  of  subroutines  that  constitute  the 
FILEGEN  processor. 

User  Response:  as  required  by  the  FILEGEN  processor. 

Comments t  A  subarea  or  soils  boundary  coordinate  file  and  a  write- enabled 
output  file  must  be  accessible  to  the  FILEGEN  processor.  The  output  file 
requires  a  logical  record  length  of  80. 
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2.  Command:  MAP GEN 

Purpose :  for  generating  subarea  or  soils  maps. 

Program  Response:  calls  the  MAPGEN  subroutine. 

User  Response:  as  required  by  the  MAPGEN  processor. 

Comments :  a  FILEGEN  produced  file  must  be  accessible. 

3.  Command :  CNCALC 

Purpose :  for  calculating  runoff  curve  numbers  and  associated  data. 

Program  Response:  call  the  CNCALC  subroutine. 

User  Response:  as  required  by  the  CNCALC  processor. 

Comments :  A  land  cover  file,  FILEGEN  produced  subarea  file,  and  either  a 

FILEGEN  produced  soils  file  or  a  soil  percentage  file,  must  be  accessible. 

4.  Command :  TTY 

Purpose :  puts  the  CNPROG  execution  in  a  "teletype"  mode. 

Program  Response:  sets  the  "TTY"  flag  to  "on". 


User  Response:  None. 
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Comments:  in  the  "teletype"  mode,  the  program  will  prompt  the  user  for 

all  input  commands  and  data. 

4.  Command :  STOP. 

Purpose :  to  stop  program  execution. 

Program  Response:  immediately  stops  program  execution. 

User  Response:  none. 

Comments :  none . 

FILEGEN  Commands 

The  FILEGEN  processor  consists  of  two  phases.  The  first  phase  is  the 
processing  of  the  boundary  point  data.  The  second  phase  divides  the  processed 
boundary  data  into  map  strips,  and  then  writes  the  scan  line  record  file. 

There  is  a  separate  set  of  commands  for  each  phase. 

Phase  I  Commands 


1 .  Command :  H 


Purpose :  this  is  a  "help"  command 
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Program  Responses  lists  the  nine  FILEGEN  phase  I  cossaands. 


User  Response;  enter  a  new  phase  I  command. 


Cowmen ts %  none . 


2 .  Command  %  STOP 


Purpose i  for  exiting  the  FILEGEN  processor  before  actual  boundary  file 
file  processing  begins. 


Program  Response:  immediately  returns  control  to  the  main  program. 


User  Responses  enter  a  neW1  MAIN  command. 


Comments ;  none. 


3 .  Command :  GO 


Purpose : 


to  perform  the  actual  processing  of  the  boundary  coordinate  file. 


Program  Response:  reads,  processes,  and  sorts  the  boundary  coordinate 
file.  The  FILEGEN  phase  II  processor  is  then  automatically  invoked. 
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User  Response;  input  a  FILEGEN  phase  II  command  after  phase  I  processing 
is  complete. 

Comments :  all  other  desired  phase  I  comments  must  be  entered  prior  to  the 
GO  command.  If  XNFILE=5,  place  the  boundary  point  cards  immediately  after 
this  GO  command. 

4.  Command:  INFILE 

Purpose :  to  specify  the  logical  unit  number  of  the  boundary  coordinate 

file . 

Program  Response:  reads  the  logical  unit  number  in  columns  1-2. 

User  Response:  input  an  integer  value  in  columns  1-2.  The  value  must  be 
right- just if ied . 

Comments  s  if  this  command  is  not  entered,  then  the  default  unit  number 
is  3. 

5.  Command ;  SHIFT 

Purpose:  shifts  the  input  boundary  coordinate  line  and  element  values 

by  specified  amounts. 


Program  Response:  reads  the  line  and  element  shift  values.  The  shift 
values  will  be  added  to  each  input  boundary  coordinate. 
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User  Response;  input  the  number  of  lines  to  shift  in  columns  1-5  and  the 
number  of  elements  to  shift  in  columns  6-10.  The  values  must  be  right- 
justified. 

1  * 

Comments ;  positive  values  will  shift  the  line  and  element  boundary  values 
down  and  to  the  right.  Negative  values  shift  up  and  to  the  left. 

6.  Command :  OUTFILE 

Purpose i  specifies  the  logical  unit  dumber  of  the  output  scan  line  refcord 
file. 

Program  Response:  reads  the  logid&l  unit  number  of  the  point  record  file. 

User  Response;  input  an  ititegei*  value  iri  columns  1-2.  The  value  must  be 
right- justified . 

Comments :  if  this  command  is  not  used,  then  the  default  unit  number  is  4. 

The  output  files  require  a  logical  record  length  of  80. 


7 .  Command :  SOILS 


Purpose ; 


specifies  that  the  boundary  coordinate  file  is  a  soils  file. 


Program  Response;  reads  the  background  soil  type  and  sets  the  flag  that 

o 

specifies  that  soils  data  is  being  processed. 
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User  Response;  enter  the  background  soil  type  (A,  B,  C,  or  D)  in  column  1. 

Comments :  if  this  command  is  not  used,  the  input  boundary  file  is  assumed 

to  be  a  subarea  file.  The  background  soil  type  is  the  type  assigned  to 
all  pixels  that  are  not  within  any  file.  The  soils  boundary  file  needs  to 
contain  the  boundaries  for  only  those  soil  types  that  differ  from  the 
background  soil  type. 

8.  Command:  DIGITIZER 

Purpose :  used  when  the  boundary  coordinates  are  digitizer  x  and  y  values 

instead  of  the  map  image  line  and  element  values. 

Program  Response:  reads  two  sets  of  coordinates;  these  sets  are  used  for 
transforming  the  input  x  and  y  values  to  the  line  and  element  values. 

User  Response:  enter  two  sets  of  coordinates;  each  set  consists  of  an 
element  (col.  1-5)  and  line  (col.  6-10)  and  the  corresponding  digitizer 
x  value  (col.  11-15)  and  y  value  (col.  16-20).  The  element  and  line 
values  are  integers.  The  x  and  y  values  are  real  and  should  include  the 
decimal  point. 

Comments :  if  this  option  is  not  specified,  then  actual  line  and  element 

values  are  assumed  in  the  boundary  coordinate  file. 
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9 .  Command :  FORMAT 

Purpose ;  to  input  the  format  specifications  of  the  boundary  coordinate 
file. 

Program  Response;  reads  the  boundary  coordinate  file  format. 

User  Response;  input  the  format  specifications,  including  the  parentheses 
in  columns  1-80. 

Comments ;  if  DIGITIZER  is  specified,  the  format  must  read  real  values 
(F  format).  If  DIGITIZER  is  not  specified,  the  format  must  read  integer 
.values  (I  format). 

The  FORMAT  command  must  be  specified  prior  to  the  GO  command. 

FILEGEN  Phase  II  Commands 


Phase  II  is  automatically  invoked  after  all  processing  from  a  phase  I  GO 
command  is  completed. 

1 .  Command :  LFBNDS 

Purpose ;  for  specifying  the  line  and  element  limits  in  the  land  cover 


file . 
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Program  Response:  reads  the  line,  element,  and  map  strip  data. 

User  Response:  several  lines  of  data  must  be  entered: 

Line  1:  enter  the  first  line  (col.  1-5)  and  last  line  (col.  6-10)  on  the 
land  cover  file. 

Line  2:  enter  the  number  of  map  strips  (max.  =  7)  in  the  land  cover  file 
(col.  1). 

Lines  3-N:  enter  the  first  element  for  map  strip  1,  the  first  element  for 
map  strip  2,  etc.,  with  one  line  for  each  strip.  Each  value  is  entered  in 
columns  1-5. 

Line  N+l:  enter  the  last  element  of  the  last  map  strip  in  columns  1-5. 

Comments :  all  values  are  integer  and  must  be  right-justified. 

2 .  Command :  GO 

Purpose :  writes  the  scan  line  record  file  to  the  unit  specified  by  the 
phase  I  OUTFILE  command. 

Program  Response:  divides  the  processed  boundary  point  data  into  map 
strips,  fills  in  the  values  between  boundary  points,  and  writes  the  scan 
line  record  file.  Control  is  then  returned  to  the  main  program. 
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User  Response:  after  all  processing  is  completed,  enter  a  new  MAIN 
command . 

Comments :  LFBNDS  must  be  entered  before  the  GO  command. 

3.  Command :  STOP 

Purpose ;  for  returning  to  the  main  program  without  using  the  GO  command. 

Program  Response;  immediately  returns  control  to  the  main  program. 

User  Response;  enter  a  new  MAIN  command. 

Comments ;  none 

MAPGEN  Commands 


The  MAPGEN  processor  produces  line  printer  maps  from  the  soils  or  subarea 
scan  line  files.  The  map  strip  line  and  element  values  of  the  MAPGEN  maps  will 
be  the  same  as  those  for  the  land  cover  map. 

1.  Command;  INFILE 


Purpose:  to  specify  the  logical  unit  number  of  the  FILGEN  produced  scan 

line  from  which  the  map  will  be  generated. 
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Program  Response;  reads  the  logical  unit  number  of  the  scan  line  file. 

User  Response:  input  the  logical  unit  number  in  columns  1-2. 

Comments :  if  this  command  is  not  entered,  the  default  unit  number  is  3. 

2.  Command:  PRINTFILE 

Purpose :  directs  map  output  to  a  logical  unit  other  than  6.  This  is 
especially  useful  when  the  user  wants  the  map  directed  to  a  line  printer 
instead  of  printing  the  map  on  the  interactive  terminal. 

i 

Program  Response:  reads  the  logical  unit  number  to  which  the  map  output 
will  be  written. 

User  Response:  input  the  logical  unit  number  in  columns  1-2. 

Comments :  if  this  command  is  not  entered  the  default  unit  is  6. 

3 .  Command :  GO 

Program  Response:  reads  the  soils  or  subarea  scan  line  record  file  and 
generates  the  map. 


User  Response:  after  the  map  has  been  written,  enter  a  new  MAIN  command. 
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Comments:  if  the  printfile  is  a  unit  other  than  6,  a  message  is  output  to 

unit  6  stating  that  the  map  has  been  written. 

4.  Command:  STOP 

Purpose:  for  exiting  the  MAP GEN  processor  before  a  map  is  generated. 

Program  Response:  immediately  returns  control  to  the  main  program. 

User  Response:  enter  a  new  MAIN  command. 

Comments:  none 


CNCALC  Commands 


The  CNCALC  processor  computes  runoff  curve  numbers  and  associated  data. 
CNCALC  requires  a  land  cover  file,  a  subarea  scan  line  file,  and  either  a 
soils  scan  line  file  or  a  soils  percent  table. 

1 .  Command :  H 

Purpose :  this  is  a  "help"  command. 


Program  Response:  prints  a  list  of  all  CNCALC  commands. 
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User  Response:  none 

Comments :  none 

2 .  Command :  FORMAT 

Purpose :  specifies  the  format  of  the  land  cover  file. 

Program  Response:  reads  the  land  cover  file  format. 

User  Response:  input  the  land  cover  file  format  (including  parentheses)  in 
columns  1-80. 

Comments :  this  command  must  be  entered  prior  to  the  GO  command.  The 
format  statement  must  be  set  up  so  that  the  line  number  is  read  as  an 
integer  followed  by  reading  a  string  of  alphameric  characters.  Each 
character  is  read  in  an  A1  format.  The  following  example  illustrates  the 
format  requirements: 

(17,  2x,  120A1 ) 

3.  Command :  SOILFXLE 

Purpose :  reads  the  logical  unit  number  of  the  soils  scan  line  record  file. 

Also,  this  command  indicates  that  CN's  will  be  based  on  a  soils  scan  line 
file  rather  than  a  soils  percent  table. 
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Program  Response:  reads  soils  scan  line  unit  and  sets  flag. 

User  Response:  input  the  soils  scan  line  unit  in  columns  1-2. 

Commands:  if  this  command  is  not  entered,  the  program  assumes  the  CN 

calculations  will  be  based  on  a  soils  percent  table. 

4.  Command :  SOILTABLE 

Purpose :  specifies  the  logical  unit  number  of  the  soils  percent  table. 

Program  Response:  reads  the  soils  percent  table  unit  number  and  then  reads 

# 

the  soils  percent  table.  Each  line  of  the  table  consists  of  the  subarea 
number  (col.  1-5),  the  percent  of  soil  types  A  (col.  6-10),  B  (col.  11-15), 
C  (col.  16-20) ^  and  D  (col.  21-25).  The  program  continues  to  read  table 
values  until  it  encounters  either  a  blank  line  or  an  end-of-file.  The 
subarea  numbers  do  not  have  to  be  in  any  particular  order.  The  subarea 
numbers  are  integer  and  the  percentages  are  real. 

User  Response:  enter  the  soils  percent  table  unit  number  in  columns  1-2. 

If  unit  5  is  specified,  the  soils  table  should  be  entered  starting  on  the 
next  line.  A  blank  line  signifies  the  end  of  the  table. 

5.  Command:  CNTABLE 


Purpose :  To  specify  the  CN  values  for  each  soil  type  and  alphameric  land 
use  symbol. 
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Program  Response:  reads  the  CN  table  on  unit  5;  each  line  consists  of  an 
alphameric  symbol  in  col.  1,  and  the  corresponding  CN  values  for  soil  types 
A  (col.  2-5),  B  (col.  6-10),  C  (col.  11-15),  and  D  (col.  16-20).  A  blank 
line  signifies  the  end  of  the  CN  table.  The  CN  values  are  real,  and  should 
include  the  decimal  point.  The  alphameric  symbols  do  not  have  to  be  input 
in  any  particular  order. 

User  Response:  enter  the  CN  table  with  the  last  line  being  blank. 

Comments :  If  this  command  is  not  entered,  then  the  default  table  shown 

in  Table  3  is  used. 

6.  Command:  LANDFILE 


Purpose ;  specifies  the  logical  unit  number  of  the  land  cover  file. 


Program  Response:  reads  the  land  cover  file  unit  number. 


User  Response:  enter  the  land  cover  file  unit  number  in  columns  1-2. 


Comments :  if  this  command  is  not  entered,  the  default  unit  number  is  4. 


7.  Command:  AREAFILE 


Purpose :  specifies  the  logical  unit  number  of  the  subarea  scan  line  file. 
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Program  Response:  reads  the  subarea  scan  line  file  unit  number. 

User  Response:  enter  the  scan  line  file  unit  number  in  columns  1-2. 

Comments :  if  this  command  is  not  entered,  the  default  unit  number  is  2. 

8.  Command :  GO 

Purpose ;  to  perform  the  CN  computations. 

Program  Response;  reads  the  scan  line  file,  computes  CN's,  and  prints 
summary  tables. 

User  Response:  after  all  processing  is  complete,  enter  a  new  MAIN  command. 

Comments ;  the  commands  FORMAT,  CNTABLE,  and  either  SOILFILE  or  SOILTABLE 
must  be  entered  prior  to  the  GO  command. 

9 .  Command :  STOP 

Purpose ;  to  exit  the  CNCALC  processor  before  actual  CN  calculations  are 
performed . 

Program  Response;  immediately  returns  control  to  the  main  program. 
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Program  Response:  reads  the  CN  table  on  unit  5;  each  line  consists  of  an 
alphameric  symbol  in  col.  1,  and  the  corresponding  CN  values  for  soil  types 
A  (col.  2-5),  B  (col.  6-10),  C  (col.  11-15),  and  D  (col.  16-20).  A  blank 
line  signifies  the  end  of  the  CN  table.  The  CN  values  are  real,  and  should 
include  the  decimal  point.  The  alphameric  symbols  do  not  have  to  be  input 
in  any  particular  order. 

User  Response;  enter  the  CN  table  with  the  last  line  being  blank. 

Comments :  If  this  command  is  not  entered,  then  the  default  table  shown 

in  Table  3  is  used. 

6.  Command:  LANDFILE 


Purpose :  specifies  the  logical  unit  number  of  the  land  cover  file. 


Program  Response:  reads  the  land  cover  file  unit  number. 


User  Response:  enter  the  land  cover  file  unit  number  in  columns  1-2. 


Comments :  if  this  command  is  not  entered,  the  default  unit  number  is  4. 


7.  Command:  AREAFILE 


Purpose:  specifies  the  logical  unit  number  of  the  subarea  scan  line  file. 
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Program  Response;  reads  the  subarea  scan  line  file  unit  number . 

User  Response;  enter  the  scan  line  file  unit  number  in  columns  1-2. 

Comments:  if  this  command  is  not  entered,  the  default  unit  number  is  2. 

8.  Command :  GO 

Purpose :  to  perform  the  CN  computations. 

Program  Response:  reads  the  scan  line  file,  computes  CN's,  and  prints 
summary  tables. 

User  Response:  after  all  processing  is  complete,  enter  a  new  MAIN  command. 

Comments :  the  commands  FORMAT,  CNTABLE,  and  either  SOILFILE  or  SOILTABLE 

must  be  entered  prior  to  the  GO  command. 

9 .  Command :  STOP 

Purpose :  to  exit  the  CNCALC  processor  before  actual  CN  calculations  are 
performed. 


Program  Response:  immediately  returns  control  to  the  main  program. 
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User  Response;  enter  a  new  MAIN  command . 


Comments:  none 


The  complete  set  of  CNPROG  commands  is  summarized  in  Fig.  15.  There  are 
a  total  of  30  commands. 


Start 


FILEGEN 


MAP GEN  CNCALC  STOP 


Phase  I: 

H 

STOP 

GO 

INFILE 

SHIFT 

OUTFILE 

SOILS 

DIGITIZER 

FORMAT 


INFILE 

PRINTFILE 

GO 

STOP 


H 

FORMAT 

SOILFILE 

SOILTABLE 

CNTABLE 

LANDFILE 

AREAFILE 

GO 


Phase  II: 
LFBNDS 
GO 
STOP 


FIGURE  15.  CNPROG  Command  Summary 
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FIGURE  17  •  EXAMPLE  LAND  COVER  MAP 
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Soils  Data  Files 

The  soils  data  can  be  input  in  two  different  ways.  The  first  way  is  to 
use  a  FILEGEN  soils  file  in  which  each  point  is  a  particular  soil  group.  The 
second  way  is  to  input  a  soil  percent  table.  The  table  contains  the 
percentages  of  each  soil  type  for  each  subarea.  The  CN  for  a  particular  point 
in  a  particular  subarea  is  based  on  a  weighted  average  CN  for  the  land  cover 
type  for  that  subarea.  The  weighted  average  CN  is  based  on  the  soil 
percentages  and  the  CN  table  values. 

Land  Cover-CN  Table 

A  table  specifying  the  CN  for  each  land  cover  and  soil  type  is  needed  for 
computing  CN's.  CNPROG  contains  the  default  land  cover-CN  table  that  is  used 
if  a  table  is  not  specified  in  the  CNPROG  input.  This  default  table  is  shown 
in  Table  3. 

TABLE  3.  Default  Land  Cover-CN  Table 


Map  Symbol 

Land  Cover  Type 

Soil 

A 

Group 

B 

C 

D 

* 

Forest 

25 

55 

70 

77 

• 

Agriculture 

36 

60 

73 

78 

+ 

Residential 

60 

74 

83 

87 

0 

Highly  Impervious 

90 

93 

94 

94 
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In  addition  to  computing  CN's  the  CNCALC  processor  also  computes  other 
desirable  information.  For  each  subarea  and  for  the  watershed  as  a  whole, 
CNCALC  computes  soil  percentages,  land  cover  percentages,  total  number  of 
classified  pixels,  and  total  number  of  pixels. 

Boundary  Point  Files 

The  accuracy  of  the  CN  estimates  is  highly  dependent  on  how  accurately  the 
soils  and  subarea  boundary  points  correspond  to  the  points  on  the  land  cover 
file.  That  is,  the  line  and  element  values  in  all  files  should  refer  to  the 
same  points  geographically. 

The  most  accurate  procedure  for  determining  boundary  points  is  to  draw  the 
boundaries  on  the  land  cover  line  printer  map.  The  boundary  point  line  and 
element  values  can  then  be  determined  accurately  and  consistently. 

The  user  can  either  input  the  boundary  point  line  and  element  values  or 
place  a  map  on  a  digitizer  and  trace  the  boundary  points  with  the  cursor.  If 
a  digitizer  is  used,  then  CNPROG  must  transform  the  boundary  points  from  the 
digitizer  x  and  y  coordinates  to  the  land  cover  map  line  and  element 
coordinates . 

The  basic  boundary  point  input  sequence  consists  of  inputting  all  of  the 
boundary  points  for  a  subarea  or  soils  region  followed  by  the  boundary  points 
for  either  a  different  region  for  the  same  subarea  or  soil  type,  or  the 
boundary  points  for  a  new  subarea  or  soil  type  region. 

Numbering  of  Subareas.  The  user  must  assign  a  unique  number  to  each 
subarea.  The  subareas  can  be  assigned  numbers  from  1  to  999.  The  subareas  do 
not  have  to  be  numbered  in  any  particular  order.  CNPROG  assigns  an  index 
value  to  each  subarea  based  on  the  order  in  which  the  subarea  boundaries  are 
entered.  These  index  values  are  used  for  computational  purposes;  all  user 
input  should  be-in  terms  of  the  user's  subarea  numbers  rather  than  the  CNPROG 
index  values. 
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Numbering  of  Soils  Areas.  The  soil  group  regions  should  be  numbered  from 
1  to  4.  The  numbers  1,  2,  3,  and  4  correspond  to  soil  group  A,  B,  C,  and  D, 
respectively.  All  of  the  regions  for  a  particular  soil  group  should  be 
grouped  together  on  the  boundary  point  file.  All  of  the  closed  regions  for 
one  soil  group  should  be  input,  then  all  of  the  closed  regions  for  a  different 
soil  type  should  be  input,  etc. 

Background  Hydrologic  Soil  Group.  The  background  soil  group  (A,  B,  C,  or 
D)  should  be  the  predominant  soil  group  in  the  watershed.  The  user  needs  to 
input  boundary  point  data  for  only  the  regions  where  the  soil  group  is 
different  from  the  background  soil  group.  The  background  soil  group  is 
specified  by  the  User  in  the  CNPROG  input. 

File  Structure 

The  boundary  point  file  structure  depends  oh  if  the  boundary  points  are 
line  and  element  values  or  digitizer  coordinates.  A  basic  difference  is  that 
the  line  and  element  values  are  input  as  integers,  whereas  the  digitizer 
coordinates  are  input  as  real  values. 

Line  and  Element  Files.  The  user  can  specify  the  format  of  the  boundary 
file  with  some  restrictions.  The  format  specification  must  read  three  integer 
values  on  each  line;  the  first  value  is  the  line  number,  the  second  value  is 
the  element,  and  the  third  value  is  either  the  user-defined  subarea  number  or 
the  soil  group  number.  The  last  boundary  point  for  a  closed  region  must  be 
the  same  as  the  first  boundary  point  for  the  region.  The  subarea  or  soil 
number  only  needs  to  be  entered  on  the  first  boundary  point  for  a  region. 
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Digitizer  Files.  The  digitizer  file  structure  is  quite  different  from  the 
line  and  element  file  structure.  The  first  line  for  a  closed  region  should 
have  the  subarea  or  soil  number  in  columns  8-9.  The  subsequent  lines  should 
contain  digitizer  x  and  y  values.  Each  line  of  digitizer  values  should 
contain  eight  sets  of  x  and  y  values.  The  end  of  input  boundary  points  for  a 
closed  region  is  indicated  in  the  input  file  by  a  zero  x  value.  The  region  is 
then  automatically  closed  in  CNPROG. 

Soil  Group  Percent  File.  A  table  containing  the  percentages  of  the 
hydrologic  soil  groups  can  be  used  in  place  of  a  soil  group  boundary  point 
file.  Each  line  in  this  file  should  contain  the  soil  group  percentages  for 
one  subarea.  The  reguired  format  for  each  line  is  given  in  Table  4. 

TABLE  4.  Soil  Group  Percent  Input  Format 


Col. 


Description  of  Contents 


1-5  Subarea  number f  integer  format 

6-25  4  five-column  real  fields  in 

which  the  first  value  is  the 
percent  of  soil  group  A, 
the  second  is  the  percent  of 
soil  group  B,  etc. 


Example  Files 

Figures  18  and  19  present  example  line  and  element  files  for  subareas  and 
soils,  respectively.  Figure  20  presents  an  example  digitizer  subarea  file. 
The  required  format  specifications  are  included  in  each  figure. 
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127  1538  5 

1  30  154  4  5 


FORMAT: 

(315) 


FIGURE  18. 


EXAMPLE  SUBAREA  BOUNDARY  POINT  FILE 
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FIGURE  19. 


EXAMPLE  SOILS  BOUNDARY  POINT  FILE 
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CNPROG  Input  Sequence 

Figure  21  shows  a  general  sequence  of  CNPROG  commands  that  can  be  used  for 
a  complete  CNPROG  analysis.  The  complete  analysis  sequence  processes  soils 
and  subwatershed  boundary  points,  produces  line  printer  maps,  and  performs  the 
CN  analyses.  The  sequence  showed  only  one  of  many  possible  sequences  because 
the  CNPROG  commands  can  be  placed  in  almost  any  desired  order  prior  to  a  "GO" 
command. 

Figures  22,  23,  and  24  are  sample  input  coding  forms  for  the  FILEGEN, 

MAP GEN,  and  CNCALC  processors,  respectively.  These  forms  can  be  used  for 
keypunching  CNPROG  input  files. 

Example  CNPROG  Runs 

Four  example  runs  are  presented  in  this  section.  The  first  three  runs 
were  performed  interactively  using  the  University  of  Maryland  UNIVAC  1100/42 
computer.  The  user  inputs  are  in  lower  case  letters  and  the  CNPROG  outputs 
are  in  upper  case  letters.  The  fourth  example  run  was  executed  in  a  batch 
mode  on  an  IBM  370  computer. 

The  runs  are  based  on  the  simplified  example  data  contained  in  Figures  16, 
18,  and  19.  The  runs  are  designed  to  illustrate  most  of  the  features  and 
options  in  CNPROG.  The  CNPROG  prompts  and  outputs  shown  in  the  examples  are 
self-explanatory. 

Table  5  shows  the  file  numbering  system  used  in  all  runs.  This  particular 
numbering  system  is  arbitrary;  the  user  must  supply  the  appropriate  computer 
system  commands  to  associate  the  unit  numbers  with  correct  files. 

The  land  cover  classifications  for  the  examples  are  shown  in  Table  6.  The 
CN’s  associated  with  these  land  uses  are  shown  in  Table  3. 
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GENERAL  CNPROG  COMMAND  SEQUENCE 
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Figure  22.  Coding  Form  for  FILEGEN  Processor 
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Figure  24.  Coding ;Form  for  CNCALC  Processor 
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TABLE  5.  Logical  File  Unit  Numbers  for  Example  Runs 


Unit  No. 

Description 

10 

Land  Cover  File 

11 

Subarea  Boundary  Point  File 

12 

Soils  Boundary  Point  File 

13 

Subarea  Scan  Line  File 

14 

Soils  Scan  Line  File 

TABLE  6. 

Example  Land  Cover  Symbols 

Symbol 

Land  Cover 

* 

• 

Forest 

Agriculture 

+ 

Residential 

0 

Highly  Impervious 
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The  first  three  runs  are  presented  in  Appendix  A.  The  UNIVAC  execution 
statements  that  start  each  run  are  underlined  to  clearly  indicate  the  start  of 
a  CNPROG  execution.  The  three  runs  together  constitute  a  complete  CN  analysis 
sequence.  The  analysis  sequence  could  just  as  easily  have  been  completed  in 
one  program  execution;  three  runs  are  used  to  illustrate  the  versatility  of 
CNPROG. 

The  first  run  generated  the  scan  line  file  for  the  soils  boundaries.  The 
FILEGEN  processor  is  used  and  then  execution  is  terminated  by  the  STOP  command 

The  second  run  uses  MAP GEN  to  produce  a  map  based  on  the  scan  line  file 
created  in  the  first  run.  Note  that  only  three  commands  are  necessary  to 
produce  the  map:  MAPGEN,  INFILE,  and  GO. 

The  third  run  is  much  more  extensive  than  the  first  two  runs.  First,  the 
FILEGEN  processor  is  used  to  produce  the  subarea  scan  line  file.  The  MAPGEN 
processor  is  used  next  to  produce  a  subarea  map.  The  CNCALC  processor  is  then 
called  twice.  The  first  call  computes  CN's  based  on  the  soils  scan  line 
file.  The  second  call  uses  a  soil  percent  table  for  computing  CN's. 

The  fourth  run  is  the  batch  mode  equivalent  of  the  first  three  runs.  The 
entire  process  of  generating  the  FILEGEN  files,  producing  the  maps,  and 
calculating  the  CN's  is  performed  in  one  execution.  The  input  file  including 
the  necessary  IBM  JCL  file  control  statements,  is  shown  in  Figure  25.  The 
resultant  output  is  shown  in  Appendix  B.  Note  that  most  of  the  CNPROG  prompts 
for  input  data  are  suppressed  in  the  batch  mode  execution. 
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//STEP1  EXEC  PGM- CMPRCG ?  REG I CM=3 0 CK 
//FTCGFCC1  EE  SYSCUT^R 

EE  BSM-SERMK C2 . CM . RREREME . ERTR ,  E I SP=SHR 
EE  ESM-SERMH C£ . CM . SC I LEME . ERTR- E I SP-SHR 
EE  ESM-SERMH  02 .  CM .  SC  I LF I LE .  ERTR  -  UM I  T-  SYS'ER? 

E I SP-  <MEU>  CRTLG>  S  BCB-  <RECFM=FEj  LRECL-S  C ,  ELKS  I  ZE=-S  1  £  O  . 
SF'RCE-  <3 120  <25  >  5>  ,  RLSE> 


//FT 1 £F 0  C 1 
//FT13FCC1 
//FT 1  OF C Cl 
// 

// 


//FT 1 1 F  C  C 1  EE  ESM-SERMH 02 . CM . RRERF I LE . ERTR , UM I T - SYSER * 

//  El  SP=^  <MEU  >  CRTLG>  ?  ECE-  <R£CFM-FB  *  LRECL^S  C  ?  ELKS  I  ZE=3 120  > 
//  SPfiCE- <3 120  <25 , 5>  >  RLSE> 

//FT 14FCC1  EE  ESM-SERMH C£ . CM . LRMETEST . ERTR ?  E I SP-SHR 
//FTC5FCC1  EE  * 

FI LEGEM 


IMFILE 

13 

CUTFILE 

1C 

FCRMRT 

<310 

SCILS 

E 

GC 

LFEMES 

ICC  125 

£ 

15CC 
i  cr^,r 
1570 
GC 

FI  LEGEM 


IMFILE 

12 

CUTFILE 

11 

FCRMRT 
<2 1 5-> 


GC 

GC 

MRPGEM 

IMFILE 

1C 

GC 

MRPGEM 

IMFILE 

11 

GC  . 


FIGURE  25 


EXAMPLE  CNPROG  INPUT  FILE 


CKCftLC 
FCRMRT 
<I5?£X? 5CR1> 
SCILFILE 
1C 

LfiNEFILE 

14 

RRERFILE  - 
11 

CMTRELE 


4 

£5 

ETicr 

70 

77 

ffl 

36 

6C 

73 

78 

c 

o  f5 

y  V 

oo 

94 

94 

4 

6C 

74 

S3 

87 

bC 

STOP 


FIGURE  25 


(Continued) 
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SUMMARY 

Previous  studies  have  shown  the  potential  of  using  remotely  sensed  data  in 
hydrologic  modeling.  Several  of  these  studies  have  focused  on  the  SCS  runoff 
Curve  Number.  The  overall  conclusion  of  these  investigations  was  that 
remotely  sensed  data  is  a  cost-effective  alternative  for  determining  the  land 
cover  required  by  the  CN  procedure. 

In  this  report  we  described  a  computer  based  procedure  for  estimating  CNs 
using  remotely  sensed  data.  The  procedure  utilizes  a  number  of  existing 
analysis  routines  and  new  procedures.  These  have  been  streamlined  for  this 
particular  application.  Descriptions  and  examples  of  each  step  in  the 
analysis  have  been  included. 
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APPENDIX  A 


OUTPUT  FOR  EXAMPLE  RUNS  1,  2,  and  3 

§ASG#A  SOILBND. 

READY  '  , 

tASG*UP  SOILFILE. 

READY 

®USE  13*  SOILBND.  .* 

EEADY 

§USE  1 OsSOILFILE. 

READY  ■  •  . 

@XQT  LIB. CNPROG 


*  CNPROG  *  AN  EASY-TO-USE  SYSTEM  FOR  COMPUTING 

*  CNPROG  *  RUNOFF  CURVE  NUMBERS 

*  CNPROG  *  CONVERTS  SPATIAL  SOILS  AND  WATERSHED  BOUNDS  TO 

*  CNPROG  *  DIGITAL  FORM  FOR  USE  WITH  DIGITAL  LAND  USE  DATA 

ENTER  ONE  OF  THE  FOLLOWING  OPTIONS: 

FILEGEN  MAPGEN  CNCALC  TTY  STOP 

TTY  .  ..  • ' '  . 

TTY 

*  CNPROG  *  TELETYPE  OPTION  SELECTED 

ENTER  ONE  OF  THE  FOLLOWING  OPT  IONS  t 
FILEGEN  MAPGEN  CNCALC  TTY  STOP 

FILEGEN  , 

FILEGEN 

*  CNPROG  *.  BOUNDARY  FILE  PROCESSOR: 

ENTER  H  FOR  LIST  OF  OPtXONS 


ENTER  BNDFIL  OPTION: 

H 

BNDFIL  OPTIONS: 

-'H  STOP  GO  INFILE  SHIFT 

OUTFILE  SOILS  DIGITIZER  -  FORMAT 

ENTER  BNDFIL  OPTION: 

INFILE 

ENTER  BOUNDARY  INPUT  FILE 

13 

UNIT  IN  COLS 

ENTER  BNDFIL  OPTION: 

OUTF  ILE 

ENTER  SCAN  LINE  FILE  OUTPUT  UNIT  IN  COLS.  1-2 

10 

ENTER  BNDFIL  OPTION: 

FORMAT 

ENTER  BOUNDFIL  FORMAT 
<3  15) 


■.ENTER  BNDF I L  OPTIONS  T  '  '  .  ■  ' /' 

soils  -  '  . 

ENTER  BACKGROUND  SOIL  TYPE  CA.»B#C*OR  D>  IN  COL®  l 

D  •  -  .  » 

ENTER  BNDF IL  OPTIONS 

GO  '  ■■■■'. 

INPUT  BOUNDARY  POINT  FILE  *  UNIT  13 
OUTPUT  SCAN  LINE  FILE®  UNIT  10 

BOUNDARY  FILE  FORMAT- <3 15 ) 

BACKGROUND  SOIL  TYPE®  D 

’INDEX  SUBAREA  NUMBER 
I  1 

1  \  I  v>: 


SOIL  POINT  FILE  CONTAINS  104  RECORDS  ■  *  •. 

.  -  ,  1  •  • 

*  CNPROG  *  NUMBER  OF  RECORDS  TO  SORT®  104 

SORT  COMPLETED  -  >  '  ’ 

.  •  *  CNPROG  *  FILEGEN  OPTIONS:  -  '■>  : 

LFBNDS  GO  STOP 

'  .  -  .  .  *  i 

ENTER  FILEGEN  OPTIONS  ‘ 

LFBNDS  ■  '  -  - .  •  .  :  ‘  ..  •  '  : 

ENTER  LAND  COVER  FILE  LINES! (COL.  1-5)  AND  L INEND (COL.  6-10) 
100  135 

\  V  •  '  •  ‘  ..  . 

' '  -  -  '  ’  "  *  'I 

ENTER  NUMBER  OF  STRIPS  IN  LAND  COVER  FILE(COL.  1)  v 

2 

ENTER  FIRST  ELEMENT  FOR  STRIP  1 (COL. I -5): 

1500  •  v  •  • 

ENTER  FIRST  ELEMENT  FOR  STRIP  2 (COL . 1 -5 ) : 

1530 

ENTER  LAST  ELEMENT  FOR  STRIP  2  (COL. 1-5)  l 
/  1570 

ENTER  FILEGEN  OPTIONS  •  /  . 

GO  ,  ■  .. 


i 


OUTPUT  FILE  UNIT®  10 

LAMB  COVER  FILE  BATA  § 

NSTRIP®  2 
LINESfa  100'- 
L2NENO*  13$ 

STRIP  1ST  ELEMENT  LAST  ELEMENT 
l  •  1500  j 529 

'  2  ...  ;  .  1530  1  §7 0 

*  CNPROG  *  SCAN  LIME  'FILE  IS  WRITTEN 
ENTER  ONE  OF  THE  FOLLOWING  OPTIONS? 


FILEGEN  . MAPGEM  CNCALG 

TTY 

STOP 

STOP 

.  STOP 

1 ' 

NORMAL  EXIT ®  EXECUTION  TIME? 

1077  “MILLISECONDS 

@XQT  LXB« CNPROG 


*  CNPROG  »  AN  EASY-TO-USE  SYSTEM  FOR  COMPUTING  • 

*  CNPROG  *  RUNOFF  CURVE  NUMBERS  • 

*  CNPROG  *  CONVERTS  SPATIAL  SOILS  AND  WATERSHED  BOUNDS  TO 

*  CNPROG  *  DIGITAL  FORM  FOR  USE  WITH  DIGITAL  LAND  USE  DATA 

'  ENTER  ONE  OF  THE  FOLLOWING  OPT  IONS  1  ,  ' 

FILEGEN  MAPGEN  CNOALG  TTY  STOP 

TTY 
TTY 

*  CNPROG  *  TELETYPE  OPTION  SELECTED 

.  ENTER  ONE  OF  THE  FOLLOWING  OPT  IONS  I  *  . 

FILEGEN  .  MAPGEN  ONCALC  TTY  STOP 

MPGEN 
MAPGEN  ' 

*  CNPROG  *  ’  MAPGEN  OPTION  ... 

ENTER  ONE  OF-  THE  FOLLOWING  OPTIONS 

INFILE  PR1NTFILE  GO  STOP 

INFILE 

ENTER  INPUT  FILE  UNIT  IN  COL.  1  **2 s 

10  v 

*  CNPROG  *  MAPGEN  OPTION  '  ' 

ENTER  ONE  OF  THE  FOLLOWING  OPTIONS 

INFILE'  PRINTFILE  Gd  STOP 

GO  •  •  ,  '  •"  '■  •"  • 


INPUT  FILE  UNIT®  10 
PRINT  FILE  .UNIT®  6 


SOILS  FILE  CONTAINS  2  STRIPS 

1ST  LINE®  100  LAST  LINE®  135 

BACKGROUND  SOIL  TYPE  IS  D 


SOILS  MAP  STRIP  I 


1500X159512510 1 15  25  XI 520  21 525 1 1 53  0  I 
1  .1  I  X  .  11-2 


too 

1 

AAAAAAAAAAAAAAAAAAAAAA 

101 

1 

AAAAAAAAAAAAAAAAAAAA 

'i  0  a 

I 

AAAAAAAAAAAAAAAAA 

S  03 

I 

AAAAAAAAAAAAAAA 

204 

I 

'•  *  .  AAAAAAAAAAAAA 

1  05 

x  -  . 

,  AAAAAAAAAAA 

106 

1 

AAAAAAAA 

S  07 

I' 

AAAAAA 

108 

1 

AAAA 

!09 

2 

AA 

110 

1 

It  1 

1 

112 

1 

113 

1 

A 

l  14 

1 

AAAA 

!!§ 

X 

-  ■  AAAAAAA 

2  16 

I  ' 

■  AAAAAAAAAA 

11? 

1 

AAAAAA AAA 

118 

I 

AAAA  AAAA  A. 

I  I  9 

I  •  ' 

’  AAAAA 

120 

1 

A. 

tat 

1 

122 

1 

123 

1 

124 

X 

•1 

125 

1 

126 

I 

127 

l 

128 

l 

BB 

229 

I  ■ 

BBB 

130 

I- 

•  ‘  .  -  BBBBB 

231 

I  ■ 

BBBBBBB 

2.32 

1 

BBBBBBBB 

233 

1 

BBBBBB3BBB 

134 

l 

BBBBBBBBBBB 

135 

1 

BBBBBBBBBBBB 

SOILS  MAP  STRIP  2 


1530 II 535 1 1540 I154S II 55 01 1555 1 1S60 1 1565  1 157 01 
l  :  i  i  I  ■  I  . ,  l  I  I  i 

I AAA AAAAAAAA AAAA AAA AAAAAAAAAAAAAAAAAAA AAAA 
XAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

..  ..  .  .  4AAAAAAAAAAAAAA 

IAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
XAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
X  AAA AA AAAA AA AA AAAAAAAA AAAA AAAA A AAA AA A A A A AA 
XAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
X  A AAA AAA AA AAAAAAAA AAAA A AAA AAA AAA AAAAAAAA A A 
XAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
X  A AAA AAAAAAAAAAA AAAAAAAA AAAAAAAA AAAAAAAA AA 
I  AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


100 

101 

102 

103 

104 

105 
S  06 

107 

108 
109 

no 
i  1 1 
1  12 
113 
1  14 

ns 
116 
S  17 
118 
I  1  9 
120 
121 
122 
123 
I  24 

125 

126 
12? 
128 
129 
S  30 

131 

132 

133 

134 

135 


AAAA AAAAAAAAAAAAAAAAAAA AAA AAAAAAAAAAAA 

'  5 ,  aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

2  AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

I  AAA AA AAAAAAAA AAAAAAAAAAAAAA AAAA 

I  aaaaaaaaaaaaaaaaaaaAaaaaaaaaa 

X  AAA AAAAAAAA AAAAAAAA AAAAAAAA 

1  AAAAA AAAAAAAAAAAAAAAAAAA 

2 A  B.'  AAAAAAAAAAAAAAAAAAA 

I  -8BB  AAA  AAAAAAAAAAAA 


.BBB- 

I  .  •'  BBBBB 

1  .  •  BBBBBBB 

1  BB3BBBBBB 

2  .  BBBBBBBBBB 

1  BBBBBBBBBBBBB 

X  BBBBBBBBBBBBBBB 
1  BBBBBBBBBBBBBBBB 
XBBBBBBBBBBBBBBBBBB 
SBBBBBBBBBBBBBBBBBBB 
2BBBBBBBBBBBBBBBBBBB 
XBBBBBBBBBBBBBBBBBBB. 
XBBBBBBBBBBBBBBBBBBBB 

IBBBBBBBBBBBBBBBBBBBBBBB 

XBBBBBBBBBBBBBBBBBBBBBBBB 

IBBBBBBBBBBBBBBBBBBBBBBBBB 


MAAAAAAAAAAAA 
AAAAAAAAAAA 
AAA AAAA 
AAAA 
A 


ENTER  ONE 
FILEGEN 
STOP 
STOP 


OF  THE 
MAPGEN 


FOLLOWING 

CNCALC 


OPT IONS l 
■  TTY  . 


STOP 


NORMAL  EXIT  ©  EXECUTION'  TIMES 


1329  MILLISECONDS 


®ASG#A  SO ILF 1LE© 
READY 

«ASG#A  AREABHDo' 
READY 

@ASG# A  LANDPILEo 

READY 

§A5G#UP  AREAFILE® 
READY 

§USE  1 0#SOILFILE® 

READY 

§USE  1 2# AREABNB® 
READY 

fUSE  14# LANDPILEo 
READY 

&JSE  1  1  #  AREAF  2LE® 
READY 

eXQT  LIB. CNPROG 


♦CNPROG  * 

*  CNPROG  * 

*  CNPROG  * 

*  CNPROG  * 


AM  EASY-TO-USE  SYSTEM  FOR  COMPUTING 
RUNOFF  GURUS  NUMBERS 

CONVERTS  SPATIAL  SOILS, AND  WATERSHED  BOUNDS  TO 

.iSIH  DIGITAL  LAND  USE  DATA 


f  ENTER  ONE  OP  THE  FOLLOWING  OPT  IONS  8 
'FILEGEN  HAPGEN  CNCALC  TTY  STOP 

TTY 
TTY 

*  CNPROG  *  TELETYPE  OPTION  SELECTED 

;  ENTER  ONE  OF  THE  FOLLOW iNG  OPT  IONS : 
FILEGEN  MAPGEN-  '  CNCALC  TTY  STOP 

FILEGEN 
FILEGEN 

*  CNPROG  *  BOUNDARY  FILE  PROCESSORS 
ENTER  H  FOR  LIST  OF  .OPTIONS 


ENTER  BNDFIL  OPTIONS 

H 

BNDFIL  OPTIONS? 

H  STOP  GO  INFILE  SHIFT 

OUTFILE  SOILS  DIGITIZER  FORMAT 

ENTER  BNDFIL  OPTIONS 

INFILE 


ENTER  BOUNDARY  INPUT  FILE  UNIT  IN  COLS.  1 -2 

12 

ENTER  BNDFIL  OPTIONS 

OUTFILE 


ENTER  SCAN  LINE  FILE  OUTPUT  UNIT  IN  COLS.  1-2 
II  -  .•  ‘  ‘  . 

,  ENTER  BNDF1L  OPTIONS  /  , 

FORMAT  '  ■ . 

'  ENTER  BQUNDF1L  FORMAT 

<3  IS)  . 

ENTER  BNDFIL  OPTION? 

GO 

•  INPUT  BOUNDARY  POINT  FILE  *  UNIT  12 
OUTPUT  SCAN  LINE  FILE®  UNIT  IS 

BOUNDARY 'FILE  FORMAT® C3 IS)  (  ' 

I 

'INDEX  ■  SUBAREA  NUMBER 

1  I  ■  ■■ 

2  2  •  . 


5  S  ’  ..  .  •...  . 

SUBAREA  POINT  FILE  CONTAINS  '  I  S3  RECORDS  ..  ' 

CNPROG  *  NUMBER  OF  RECORDS  TO  SORT*  183 
SORT  COMPLETED.  • 

*  CNPROO  *  .PILR6EN  OPT  IONS  » 

'  LFBNDS  .  60  STOP  '  ■  '  '  .  ■  •  .  •'  ‘  ■ 

ENTER  FILEOEN  OPTIONS  .  ■  : 

LFBNDS  -  : 

ENTER  LAND  COVER  FILE  LINEST (COL.  1-5)  AND  L INEND (COL.  6-10) 
180  135 

ENTER  NUMBER  OF  STRIPS'  IN  LAND  COVER  FXLE(COL.  1) 

2  ■  .  •  • 

ENTER  FIRST  ELEMENT  FOR  STRIP  l(C0L®l“5)s 
150  0  ■'  '  . 

ENTER  FIRST  ELEMENT  FOR  STRIP  2 (COL. 1-5): 

1530 

ENTER  LAST  ELEMENT  FOR  STRIP  2CC0Lot-5): 

1570 

ENTER  FILEGEN  OPTION?  '  ’ 

GO 


OUTPUT  FILE  UNIT®  n 

.  '  •  ’  V  •*  ••  .* 

LAND  COVER  FILE  DATA s  / 

-NSTRIP®  2 
L  ICIEST®  100 
L INEND®  135 

STRIP  '  1ST  ELEMENT  LAST  ELEMENT 


. 

.  2 

2  ; 

1500 

IS30 

*•  1529  • 

1570  ' 

'  \  ... 

CNPROG 

«- 

AREAS . 

2 

AND  , 

3 

OVERLAP 

AT 

COLUMN 

1537 

LINE 

1  05 

# 

CNPROG 

AREAS 

2 

AMD 

3 

OVERLAP 

AT 

COLUMN 

1537 

LINE  _  . 

106 

*. 

CNPROG 

AREAS 

2 

AMD 

3 

OVERLAP 

AT 

COLUMN 

1537 

LINE  ! 

107 

CNPROG 

AREAS ' 

2 

AND 

3 

OVERLAP . 

AT 

COLUMN 

1536 

LINE 

1  10 

'* 

CNPROG 

.  AREAS  ■ 

2 

AND 

'3 

OVERLAP 

AT 

COLUMN 

1536 

LINE 

111 

*  CNPROG  *  SCAN  LINE  FILE  IS  WRITTEN 


ENTER  ONE  OF  T HE  FOLLOWING  OPTIONS! 


F1LEGEN 

MAPGEN 

•  MAPGEN 

CNCALC 

. TTY  STOP 

MAPGEN 

•  '  ’ 

■  *  CNPROG 

'*  MAPGEN 

OPTION 

ENTER  ONE 

OF  THE  FOLLOWING 

OPTIONS 

INF ILE 

INFILE 

PR.IMTFILB 

-  GO 

STOP 

ENTER  INPUT  .FILE  UNIT  IN  GdL.  1-2: 

11. 

, 

*  CNPROG 

*  MAPGEN. 

OPTION 

■  ENTER  ONE 

OF'  THE  FOLLOWING 

OPTIONS 

INFILE 

PRXNTFILE 

GO 

STOP 

GO 


INPUT  FILE  UNIT®*  11 
PRINT  FILE  UNIT®  6' 


SUBAREA  FILE  CONTAINS  2  .STRIPS 
1ST  LINE®  10.0  LAST  LINE®  135 


INDEX  SUB-AREA  NO®  MAP  SYMBOL 

1  '  I  A 

2  2  B 

3  3  G 

4  *  ■  4  .  •*  D  ■ 

5  '  '  ,  S'  ’£  •  ' 

SUB- AREA  MAP  STRIP  NUMBER  1 


1500X15051 151  OX  15 15X1520 1 1525X1530  I 
X  .1.1.  I  I  XX 


100  1 
10.1  X 

102  X 

103  I 
i  04  I 
I  05  I 

106  I 

107  I 

108  1 
109  1 

no  i 
in  x 
112  1 
11.3  l 

114  X 

115  I 

116  1 
S17  1 
118  l 
ISf  1 

120  I 

121  I 

122  I 

123  I 

124  1 

125  S 

126  X 

127  1 

128  X 

129  I 

130  I 

131  I 

132  I 
1.33  1 

134  2 

135  1 


A  '  B 

A AAA  BBBBBBBB 

AAAAAAABBBBBBBBBBBBBBBB 
AAAAAAABBBBBBBBBBBBBBBB 
AAAAAAABBBBBBBBBBBBBBBB 
AAAAAAAABBBBBBBBBBBBBBBB 
•  AAAAAAAABBBBBBBBBBBBBBBB 
AAAAAAAABBBBBBBBBBBBBBBB 
AAAAAAAAABBBBBBBBBBBBBBBB 
AAAAAAAAABBBBBB88BBBBBBBB 
.  AAAAAAAAABBBBBBBBBBBBBB  B 
. ‘ AAAAAAAAABBBBBBBBBBBBBDDD 
AAAAAAAAAABBBSBBBBB  DDDDDD 
AAAAAAAAAABBBBBB  DDDBDDDBD 
AAAAAAAAAABBBB  DDDBDDDDBDD 
• AAAAAAAAAAAB  BDDDDDBOBBDBBD 
AAAAAAAAAAABBDDBBBDBDDBBDDD 
AAAAAAAAAAA  BDBBBDDDBD 

AAAAAAAAAAAAAAAAAAAA  ’  DDDB 
AAAAAAAAAAAAAAAAAAAAAAAAAAAA 
AAAAAAAAAAAAAAAAAAAAAAA 
AAAAAAAAAAAAAAAAA 
AAAAAAAAAAAA 
AAAAAA 
A 


SUB- AREA  MAP  STRIP  NUMBER  2 


.  1530 1 1535X1540  XI 545  XI 5S0  X 155511560 11565 1 1570 1 
III  X  X  1  1  •  X  I 

100  X 

1011  '  CGCGCGO 

102  I  '  GCGCGCGGGG 

103  X.  ■  GCGGGGOGCCOCC 

104  X  CCCGCGCGGGCG  ' 

105  I  CCCCCCCCCCCCCDDDDDDDDDDD 

106  IBpBBBBBCCCGGGGGGGGGGDDDDDDDDDDDD 

107  XBBBBBBBCGCCGCCGGGOGGBBDDDDDDDDBB  • 

108  IBBBBBBBGCCGGGGCGGGGOBBDBBBBBDDDDB 

109  XBBBBBBBCCCCGGGGGCCCCODDDDODQDDOOD 

110  IBBBBBBGCGGCOGCGCGGCGBDBBDDBDDBDDDD 

1 1 1  XBBBBBBGGCGOGCC  DDDDDDBDDDDDDDDDDDD 
.112  XBBBBBBC  DDDDDDDDDDDDDODDDODDDDOD 

1 1 3  IBBBBBBDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 

114  IBBB  DDDDDDDDDDDDDDDDDDDDDDDDODDDDDDD 

11 5  XBODDDDODBBDDOODDDDBDDOBDDDDDDDODDDDD 

116  IDDDDODDDOODDDDOODDDOODDDOPDOODDDDDDOD 

117  IDDDDDDDDDDDDDDDODBDDDDBDDDDDDDDDDDD  . 

118  IDDOODDDDDDDDOpODDDDObODDODODDODDDD 

119  IDDDDDDDDDDOBDDDDODDDDBBOBDBDDDDD 

120  IDDDDDBDDDDDBBDDDBDDBBBDBBDDDDD 

121  XODDDOODDDDDODDDDDOOODDDDDDDDO 

122  XDDBDBBDDBDDDBDDBBDDDBBBBDBB 

123  IDDDDDDDDDDDDDDDDDDDDDDDDD 

124  XDDDDODDBBDODBOBDBODBBBD 

125  IE  BDDBBBDDBBBBBDDBBDDD  1 

126  IEEEEE  OBBBOBDBODDB 

127  XEEEEEEEEE  BBDDBDDB 

128  XEEEEEEEEEEE  BBBDD  ■■  l  ' 

129  IEEEEEEEEEEEEE  D 

130  IEEEEEEEEEEEEEEE 

131  IEEEEEEEEEEEEE 

132  X  EEEEEEE 

133  X  E 

134  X 

135  X 

ENTER  ONE  OF  THE' FOLLOWING  OPTIONS* 

FILEGEN  MAPGBN  CMGALC  TTY-  STOP 

CNCALC 
CNCALG 

*  CNPROG  *  'CNCALC  OPTIONS 

ENTER  H  FOR  LIST  OF  CNCALG  OPTIONS 

ENTER  CNCALC  OPTIONS  , 

H 

*  CNPROG  *  LIST  OF  ONCALC  OPT  IONS  * 

FORMAT  SOILF1LE  SOXLTABLE  CNTABLE 
LANDFXLE  ARBAFILE  GO  H  STOP 


ENTER  CNOALG  .OPTIONS 

FORMAT  ' 

ENTER  FORMATS 

a5#2X#50AI3>  .  ’  •  . 

ENTER  CNGALO  OPTIONS 
SOILFILE 

..  ENTER  SOIL  FILE  UNIT  IN  COL®  1-2 

10 

ENTER  CNCALC  OPTION 8  r  '  • 

LANDP ILE  -  ' 

,  ENTER  LAND  COVER  FILE  UNIT  IN  COL*  1-2: 
14  .  .  '  ... 

ENTER  CNGALG  OPTIONS 

AREAFILE  ,  ' 

ENTER  SUB-AREA  FILE  UNIT  IN  COL.  1-2 

II  ...  ' 

-ENTER  CNCALC  OPTION! 

CNTABLE 


ENTER  TABLE  OF  VALUES  FOR  LAND  USE  SYMBOL (COL.  1  >, 

CM  FOR  SOIL  TYPES  ACCOL*  2-5),  BCCOL*  6~10>,  CCCOL*  11-15), 

AND  D (COL®  16-20 ) 

*  25  • 55  70  11 


.  36 

60  73  78 

0  90 

93  94  94 

+  60 

74  83  8? 

/  ’ 

ENTER 

CNCALC  OPTIONS 

GO 

.  CNCALC  FILE  DATA : 

a  : 

LAND  USE  FILE  UNIT®  14 

SUB-AREA  FILE  UNIT® 

SOIL 

SCAN  LIME  FILE  UNIT®  10 

LAND 

COVER  FILE  FORMAT 

® (I5,2X,50Al ) 

CN 

»  LAND  USE  -  SOIL 

TABLE: 

SYMBOL 

A  B  C 

D 

& 

25.0  55©  0  70  &  0 

77®  0  • 

36 © 0  60g0  73  ©  0 

7  8o  0 

0 

90<=0  93  e0  94  4 

f  4  o  0 

'  + 

.  60o0  74*0  83 © 0 

87  ®  0 

SUB-AREA  FILE  INDEX  VALUES: 


INDEX 

1 

2 

3 

4 

5 


SUB -AREA  NGo 

I 

•  2 
.  3 

4 

5  ' 


SUB-AREA  FILE  BOUNDS §  -  .  . 

FIRST  LINE®  .100  LAST  LIME®  135  ;  ■ 

ELEMENT  LIMITS?  •  ..  '  /  .  '  '  ,  .  -  ’ 

'  STRIP  1ST,  ELEMENT  .LAST  ELEMENT 

’  '  1  IS0§  ,  .  '  .  1529  < 

2  •  '  153©  '  1570  ’  .  ; 

BACKGROUND  SOIL' TYPE  IS  0 

"  .  •  ■  .  .  ‘ 

‘  .  •’  .  :  ...  :  ,  \  •  .  • 

'  CN  SUMMARY  TABLE t 

INDEX  '  SOB-AREA  TOT  ©  PIXELS  CLASSIFIED  PIXELS  CN 


1 

8 

256* 

256. 

'  •  80 

.•  :  2  . 

239*  ■ 

239. 

71 

3  ;'V 

3 

I30o 

130.  ' 

26 

'  4 

*.■  ■;  4  • 

607,* 

607. 

65 

5 

,  ■  \5 

75© 

'  '  ’ 

75. 

65 

TOTALS! 

\\'1 

. .  :  >  • 

( 

'  1307a' 

1307. 

•  '  65 

SUB-AREA 

-  SOIL 

PERCENT 

TABLE : 

SUB-AREA 

•A 

'  B 

G 

D 

2 

■  o  0  0 

2*34 

©00 

97.66 

2 

33.47 

'  ©00 

©00  •’ 

66.53 

,  3 

I00O00 

•‘©00 

©00 

•  00 

'  4 

.47*94 

E4©33 

©00 

37.73 

•  5 

',.00 

96©00 

'  ©00 

4*00  ' 

,  TOTALS l 

38*33 

.12*62 

©  0  0 

49*04 

LAND  COVER  PER''  T/: :  EE 


LAND  COVER  SYMBOL 


SUB-AREA 

0 

O 

4= 

1 

i00 

?6®56 

4.69 

1  8o75 

2 

30  ©54 

24*69 

l  l  ©72 

33*05 

3 

93.08 

6©  15 

o0  0- 

oil 

4 

12*36 

27  ©  02 

29®  98 

30  ®64 

5 

•-  ©00 

6S  o  33 

I  ©  33 

33  o  3  3 

TOTALS  § 

20*58 

36©42 

I  7  ©  06 

25  ©  94 

END  OF  CNCALC  ROUTINE 

'  ENTER  ONE  OF  THE  FOLLOWING  OPTIONS! 
FILEGEN  MAPS EM  CNCALC  TTY  STOP 

CNCALC 
CNCALC 

*  CNPROG  *  CNCALC  OPTION i 
ENTER  H  FOR  LIST  OF  CNCALC  OPTIONS 
ENTER  CNCALC  OPTION? 

H 


*  CMPROG  *  ■  LIST  OF  CNCALG  OPTIONS* 

FORMAT  SOILFILE  .  SQILTABLE  C NT ABLE 

LANDFILE  ARBAPSLE  GO '  . H  .  STOP 

-  ENTER. CNCALG  OPT SON 8 

FORMAT  . 

ENTER  FORMATS  ' 

CX5# 2X#50Ai )  „ 

ENTER  CNCALG  OPTIONS 
'SOILFILE 

ENTER  SOIL  FILE  TO IT  IN  COL©  1-2 

10 

ENTER  CNCALG  OPTION  § 

SQILTABLE 

ENTER  SOIL  Vile  TABLE  UNIT  IN  GOL®  i-2* 

S  •  ■  ■ 

ENTER  SUB-AREA  -  SOIL  PERCENT  TABLE  * 

1  0*0  2*34  0.8,97*66  .  •  ;  - 

233*47  8.©  0  o866©53 

3  100®  8  ®  §  8  e  0  0  ©  0  ; 

4  <  ' 

447. 94/4 ®33« 

447.9414*3  O®0  37©7« 

447 .'94  14 ©33  0.8  4©0 

5  Q  ©  0  1  96  ©  @©0  4  o0 

ENTER  CNCALG  OPT  I  ON  8  •:  '  .  - 

SQILTABLE 

ENTER  SOIL  FILE  TABLE  UNIT  IN  COL.  1-2* 

5 

ENTER  SUB-AREA  -  SOIL  PERCENT  TABLE* 

I  0*02*34  0  ©  ©97  ©66 

233*47'  0.0  0*066*53 

3100®  0.  ®  0  0o0  flo8 

447.9414*33  O©037©73 

5  8*0  96  ®  IcB  4o| 

ENTER  CNCALG  OPTIONS 
LANDFILE 

.  ENTER  LAND  C0UER  FILE  UNIT  IN  COL.  1-2* 

14  -0  '  ,  .  : 

ENTER  CNCALG  OPTIONS  '  . 

AREAFXLE  ’  ' 

•  •  ENTER  SUB-AREA  FILE  UNIT  IN  COL©  1-2 

H  -  ' 

ENTER  CNCALG  OPTIONS 
GO 


CNGALC  FILE  BATA!  '  .  •  •  - 

LAND  USE  FILE  UNIT®  14'.  SUB-AREA  FILE  UNIT®  11 
SOIL  TABLE  FILE  UNIT®  S 


LAND  COURR  FILE  FORMAT®  C  15,9 2X*5 ©A!  > 


:  ,  ON  - 

LAND  USE 

:  -  SOIL  TABLE i 

SYMBOL 

A  B 

C 

'  D 

* 

25  o  0  55  © 

0  70  a@  ■ 

lie®.  , 

o 

■36  ©  0  6  0  o 

©  73  ©  0 

IBo® 

0 

90 4  93© 

§  94o;e, 

94.0 

-  + 

S'©  ©0  7,4  • 

0  S3  o0 

87  ©0 

"  SUB-AEEA 

-  SOIL 

PERCENT 

TABLE 

SUB -AREA 

•  A  :■ 

B 

0 

D 

1 

■  o00 

,2®34 

.  ©0® 

97  ®  66 

2 

@©e 

.00 

S6©53 

3 

'  100.00 

©00 

®  0  0 

.00 

=-.  4 

4?  o  94 

14.33 

0  0  0 

37  o?3 

■  5 

.  •  ©  0  0 

96.00 

»00 

4.00 

SUB-AREA  FILE  INDEX  ^ALUES* 


INDEX 
S 

2  • 

3 

4  ' 

5 

SUB-AREA  FILE  BOUNDS  § 

FIRST  LINEs  -100  LAST  LINE-  135 

ELEMENT  LIMITS! 

STRIP  1ST  _  LAST  ELEMENT 

1  .1500  ’  •  1529 

2  1530  ’  1570 


SUB*AREA  NOo 
1  ''  ' 

'  ’..'2 

;  '3 

4' 

5 


CN  SUMMARY  TABLE 'I 
INDEX  SUB-AREA 

1  '  I 

2  2 

3  3 

4  -  4 

5  5 


( 


.  PIXELS ' 

CLASSIFIED  PIXELS 

CN 

256  © 

256. 

80. 

23-9© 

•  239. 

71  . 

130© 

130. 

26  * 

607  © 

607. 

71  . 

75  o 

75. 

66. 

1307© 

1307. 

68. 

TOTALS  s 


LAND  COVER  PERCENTAGE  TABLE  % 

LAND  COVER  SYMBOL 

■SUB-AREA  ®  ■  0 .  + 

>  I  41  76*56  4*69  1 S®  75 

2  30oS4  24®S9  !  S  ©72  33.05 

3  .  93®iS  6® IS  ®0O  oil  . 

4  .  12o36  27©02  29.98  30*64 

5  v  .00  SS®33  31 .33  .33*33 

TOTALS  %  20*58  36*42  I7©06  25*94 

END  OF  CNCALG  ROUTINE 

ENTER  ONE  OF  THE  FOLLOWING . OPT IONS s 
PILEGEN  .  MAP6BN  ONCALC  TTY 

STOP 

STOP  "  -  .  . 

NORMAL  EXIT.  EXECUTION  TIME? 


/ 


STOP 


7 S 69  MILLISECONDS. 


APPENDIX  B 


OUTPUT  FOR  EXAMPLE  RUN  4 


*  CNPROG  *  AN  EASY-TO-USE  SYSTEM  FOR  COMPUTING 

»  cnprog  *  runoff  curve  numbers 

*  CNPROG  *  CONVERTS  SPATIAL  SOILS  AND  WATERSHED  BOUNDS  TO 
»  CNPRQG  *  ■  DIGITAL  FORM  FOR  USE  WITH  DIGITAL  LAND  USE  DaTa 

enter  one  of  the  following  options* 
filegen  mapgen  cwcalc  tty  stop 

FILEGEN 

»  CNPROG  *  BOUNDARY  FILE  PROCESSOR 

INPUT  BOUNDARY  POINT  FILE  =  UNIT  13 
OUTPUT  SCAN  LINE  FILE-  UNIT  10 

BOUNDARY  FILE  FORMAT^ C 315 ) 

BACKGROUND  SOIL  TYPEs  D 

INDEX  SUBAREA  NUMBER 

1  1 

11 
2  2 

soil'  point  FILE  CONTAINS  104  Pecoros 

*  CNPROG  *  NUMBER  OF  RECORDS  to  SQRT  =  104 

SORT  COMPLETED 

»  CNPROG  *  FILEGEN  PHASE  II  PROCESSOR 


OUTPUT  FILE  UNIT-  10" 

land  COVER  file  DA  I  a  © 

NSTRIP=  2 

LI  NEST  =  100 
L  INENDa:  135 

STRIP  1ST  ELEMENT  LAST  ELEMENT 

1  1500  1529 

2  1530  1570 

*  CNPROG  *  SCAN  LINE  FILE  IS  WRITTEN 

enter  one  of  the  following  options. 

FILEGEN  .  hAPGEN  CNCALC  TTY  STOP 
FILEGEN 

*  CNPROG  *  'BOUNDARY  FILE  PROCESSOR 

INPUT  BOUNDARY  POINT  FILE  =  UNIT  12 
■OUTPUT  SCAN  LINE  F 1LES  UNIT  11 

BOUNDARY  FILE  FOR|*AT=  <315 ) 

INDEX  SUBAREA  NLM0ER 

1  1 

2  2 


5 


5 


SUBAREA  POINT  FILE  CONTAINS  183  RECORDS 

*  CNPROG  *  NUMBER  OF  RECORDS  TO  SORT=  lo3 

sort’  completed 

*  CNPROG  *  FILEGEN  PHASE  II  PROCESSOR 


OUTPUT  FILE  UNIT  =  II 


LAND  COVER  FILE  DaIAo 

NSTRIP=  2 

LINEST=  100 
LINEND=  135 


STRIP  1ST  ELEMENT  LAST  ELEMENT 

1  1500  1529 

2  1530  1570 


CNPROG 

«• 

AREAS 

2 

AND 

3 

OVERLAP 

AT 

COLUMN 

1537 

LINE; 

105 

CNPROG 

«■ 

AREAS 

2 

AND 

3 

OVERLAP 

at 

COLUMN 

1537 

LINE 

106 

# 

CNPROG 

«• 

AREAS 

2 

AND 

3 

OVERLAP 

AT 

column 

1537 

LINE 

107 

* 

CNPROG 

AREAS 

2 

AND 

3 

OVERLAP 

AT 

COLUMN 

1536 

LINE 

no 

CNPROG 

* 

AREA? 

2 

AND 

3 

OVERLAP 

AT 

COLUMN 

1536 

LINE 

111 

»  CNPROG  *  SCAN  LINE  FILE  IS  WRITTEN 


ENTER  ONE  OF  THE  FOLLOWING  OPTIONS. 
FILEGEN  MAPGEN  CNTCALC  TTY  STOP 

MaPGEN 

*  CNPROG  *  MAPGEN  PROCESSOR 


INPUT  FILE  UNIT-  10 
PRINT  FILE  UNITs  6 


SOILS  FILE  CONTA 
1ST  LI NE=  100 


INS  2  STRIPS 
LAST  LINE= 


BACKGROUND  SOIL  TYPE  IS  0 


135 


15001 15051 15101 15151 1520  1 15251  15301 
I  I  I  I  I  I  I 


100 

I 

aaaaaaaaaaaaaaaaaaaaaa 

101 

I 

aaaaaaaaaaaaaaaaaaaaaa 

102 

1 

aaaaaaaaaaaaaaaaaaaaaa 

103 

I 

aaaaaaaaaaaaaaaaaaaaaa 

104 

I 

aaaaaaaaaaaaaaaaaaaaaa 

105 

I 

aaaaaaaaaaaaaaaaaaa 

106 

I 

aaaaaaaaaaaaaaaaa 

10? 

I 

AAAAAAAAAAAAAA 

108 

I 

A AAAAAAA A AA 

109 

I 

AAAAAAAA 

110 

I 

AAAAAA 

111 

I 

AAA 

112 

113 

I 

I 

AAAAAAAA 

114 

I 

AAAAAAA AA 

115 

I 

AAAAAAAAAA 

116 

I 

A A AAA A AAA A 

117 

I 

AAAAAAAA 

118 

I 

AAAAAA 

119 

I 

AAA 

120 

I 

A 

121 

122 

123 

124 

125 

126 

127 

128 

I 

I 

I 

I 

I 

I 

I 

I 

BB 

129 

I 

BBS 

130 

I 

8BBBB 

131 

I 

BBBBBBB 

132 

I 

BBBBBBBB 

133 

I 

BBBBBBBBBB 

134 

I 

BB8BBBBBBBB 

135 

I 

BdBBBBBBBBBB 

153011535115401154511550115551156011565115701 
I  I  1  I  I  I  I  I  I 
1 AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
I AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
I AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
I aaaaaaaaaaaa AA AAAAAAAAAAAAAAAAAAAAAmAA aaa 
I AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAaAAAAA 
IAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAmAAAAA 
IAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
I AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAMAAAAA 
IAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
I AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
I AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
IAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
lAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
I A  AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
AAAA AAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
AAAAAAAAAAAAAAAAAAAAAAAAAaAAAAa 
AAAAAAAAAAAAAAAAAAAAAAAAAAA 
baaaaaaaaaaaaaaaaaaamaaa 


100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 


B8B 
bbbbb 

8BBB8BB 
BBBdB80B8 
BBBBeBBBBB 

bbbbbbbbbbbbb 
bbbbbbbbbbbbbbb 
bbbbbbbbbbbbbbsb 
ibbbbbbbbbbbbbbbbbb 

IBBBBBBBBBBBcBBBBBBB 

ibbbbbbbbbbbbbbbbbbb 

iBBBBBBBBBBdfcBBBBBBB 

IBBBBBBBBBBBcBBBBBBBo 

ibbbbbbbbbbbbbbbbbbbbb 

ieBBBBBBBBBBfcBBBBBBBbBBb 
iBBBBBBBBBBBfcBBBBBBBBBBBB 
I  BBBBBBBBBBBcB BBBBB Bb BBBBB 


AAAAAAAAAAAAAAAAAAA 

AAAAAAAAAAAAAAaA 

aaaaaaaaaaaa 

aaaaaaaa 

aaaaa 


enter  one 
filegen 
mapgen 


OF  THE 
mapgen 


following 

cncalc 


options, 

TTY 


STOP 


*  CNPROG  *  MAPGEN  PROCESSOR 


input  file 

PRINT  FILE 


UNIT  = 
UNI  T - 


11 

0 


SUBAREA  FILE  CONTAINS 

2  STRIPS 

1ST  L I NES  100  LAST 

L I Nc=  135 

index  SUB-AREA  no* 

MAP  SYMBOL 

1  1 

A 

2  2- 

B 

3  3 

C 

SUB^AREA  MAP  STRIP  NUMBER  1 


15001150511510115151152011525115301 

I  I  I  I  I  I  i 


100  I 

101  I 

102  I 

103  I 

104  I 

105  I 

106  I 

107  I 

108  I 

109  I 

110  I 

111  I 

112  I 

113  I 

114  I 

115  I 

116  I 

117  I 

118  I 

119  I 

120  I 

121  I 

122  I 

123  I 

124  I 

125  I 

126  I 

127  I 

128  I 

129  I 

130  I 

131  1 

132  I 

133  I 

134  I 

135  I 


A  B 

AAAA  B6BBBBBB 

AAAAAAABBBBddBBBBBBBBBB 
AAAAAAABBBBddBBBBBBBBBB 
AAAAAAABBBBddBBBBBBBBBB 
AAA AAA AABBBBdtiBBBBBBBBBB 
AAAAAAAABBBBdBBBBBBBBBBB 
AAAAAAAABBBBddBBBBBBBBBB 
A AAAAAAAABBBBdBBBBBBBBBBB 
AAAAAAAAABBBBBdBriBdBBBBBB 
AAAAAaA AA0006ddBB3BBBBB  D 
AAAAAAAAABBBBBdBBBBBBDDDD 
AAAAAAAAAABBBSddBdB  DDDDDD 

aaaaaaaaaabbbbbb  dqddddddd 
aaaaaaaa aabbbb  dddoodddddd 
aaaaaaaaa aab  dooodddodddodd 
aaaaaaaaa a adddooodddddddddd 
aaaaaaaaa aa  dodododddo 
aaaaaaaaaaaaaaaaaaaa  ddod 

AAA AAA AAAA A AAAAAAA AAAAAAAA A A 
AAAAAAAAAAAAAAAAAAAAAAA 
AAAAAAAAAAAAAAAAA 
A A AAAA AAAA A A 
AAAAAA 
A 


15301153511540115451155011555115601156511570 
I  I  I  I  I  I  I  1 


100  I 

101  I 

102  I 

103  I 


CCCCCCC 

cccccccccc 

ccccccccccccc 


104  I  cccccccccccc 

105  I  cccccccccccccdododdddddo 

106  IB6B8BBBCCCCCCCCCCCCCD0DDDDD00D00 


I 

I 


107  IB6B8BBBCCCCCCCCCCCCC00DDDDDD0D00 

108  IB86B680CCCCCCCCCCCCCDDDDDDDDDDODD 

109  ibbbbbbbcccccccccccccoododdddddddd 


110  IBBBBBBCCCCCCCCCCCCCCOOQDDDDDDDDDDD 

111  IBBBBBBCCCCCCCC  DDODODODODDDDDDDODO 

112  ibbbbbbc  DODDODODDDODDDDODDODODDO 

113  IBBBBBBDODDOCDODDDODOODDODDDDDDOODDO 


114  I BBB  DODODDDCDDDDQQDUDDDODDDDDDDDDDDO 

115  IBDODDDDDDDOCDDDDDDODODDDDDDDOOODDDDU 

116  IDODDDDDODODCDDODDDDODDDODOODDDDODDOUD 

117  IDODDDDDODOOCODDDOODDDDDDDODDDOODDDIJ 

118  IDDDDDODODOOODDDDODQOOOODDDDODDOODO 

119  IDDDDDDDDDDDCDDDDODDODDDDDDDDDDDQ 

120  idddddddddddcdddddddodddododddd 

121  iddddddddddddddddqddoddddddddd 

122  iddddddddddodddddddoddddoddd 

123  ioddddddddddcddodddddddddd 

124  lOODODOOODDOCDDDOOOubODD 

125  IE  ODOODODObODDODDOOODD 

126  IEEEEE  DDDCDDDDbOOD 
127.  IEEEEEEEEE  OCDOODDD 

128  IEEEEEEEEEEE  DDDob 

129  IEEEEEEEEEEEEE  0 

130  IEEEEEEEEEEEEEEE 

131  IEEEEEEEEEEEEE 

132  I  EEEEEEE 

133  I  E 


134  I 

135  I 


enter  one  of  the  following  options. 

FILEGEN  MAP6EN  CNCALC  TTY  STOP 
CNCALC 


'*  CNPROG  *  CNCALC  PROCESSOR 

FORMAT 

soilfile 

landfile 

areafile 

CNTABLE 

go';; 

CNCALC  FILE  DATA. 
sub-area  FILE  UNIT=  11 

LAND  USE  FILE  UNITS-  14  14  14  14  14  14  14 
SOIL  SCAN  LINE  FZLc  UNIT  =  10 

LAND  COVER  FILE  F CRMAT=  ( 1 5 » 2X ♦ BOA  1 ) 


CN 

-  LAND 

USE  - 

SOIL 

TABLE 

SYMBOL 

A 

B 

C 

0 

* 

25  s  0 

55  o  0 

?0.0 

77.0 

36.0 

60.0 

7  3.0 

78  o  0 

/a  s\  rv  *n  ft  L  /.  A  iT  /.  . 


INDEX 

1 

2 

3 

4 

5 


SU8-AREA  N0. 
1 

2 

3 

4 

5 


SU8"AR£A  FILE  BOUNDS. 

FIRST  LIN£s  100  LAST  LXNE=  135 


element  limits. 

STRIP  1ST  ELEMENT  LAST  ELEMENT 


1  -1500  1529 

2  1S30  1570 


background  soil  type  is  d 

CN  SUMMARY  TABLE  . 

index  sub-area  tot.  pixels  classified  pixels  cn 


1 

i 

256. 

256. 

80. 

2 

2 

239. 

239. 

62. 

3 

3 

130. 

130. 

26. 

4 

4 

607. 

607. 

64. 

5 

5 

75. 

75. 

65. 

TOTALS 

1307. 

1307. 

63. 

SUB-AREA 

-  SOIL 

PERCENT 

TABLE. 

SUB-AREA 

A 

B 

C 

D 

1 

-  0.0 

2.34 

0.0 

97.66 

2 

55 » 65 

0.0 

0.0 

44.35 

3 

100.00 

0.0 

0.0 

0.0 

4 

SO. 91 

14.33 

0.0 

34.76 

5 

0.0 

96.00 

0.0 

4.00 

TOTALS. 

43.76 

12.62 

0.0 

43.61 

land  cover  percentage  table. 


LAND 

i  COVER 

symbol 

AREA 

» 

• 

0 

♦ 

1 

0.0 

76.56 

4.69 

18.75 

2 

30  ®  54 

24.69 

11.72 

33.05 

3 

93.08 

6.15 

0.0 

0.77 

4 

12.36 

27.02 

28.17 

32.45 

5 

0.0 

65.33 

1.33 

33.33 

ggwa 


TOTALS.  20»5d  36»<i2  16 »cd  cb. 
END  OF  CNCALC  ROUTINE 


cf.NJr?M°Ne  2au-cm  FCLL0wIN6  OPTIONS. 
FILEGEN  MARTEN  CNCALC  TTY  STOP 

STOP 


I 


I 
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